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MATHEMATICAL BIOPHYSICS OF CELL RESPIRATION II 


H. D. LANDAHL 
The University of Chicago 


The recent approximate method developed by N. Rashevsky in 
his book Mathematical Biophysics is shown to lead to results equiva- 
lent to those obtained in the preceding paper on cell respiration. 
An estimate is made of the value of the total diffusion resistance 
to lactic acid for unfertilized Arbacia eggs. Several more general 
reactions are treated to illustrate the effect of intermediate products. 
Expressions for the glycolytic coefficient are derived in terms of the 
parameters of the cell. 


I 


In continuing the theoretical study of cells with particular ref- 
erence to respiration, we shall first show that the method recently 
outlined by N. Rashevsky (Rashevsky, 1938) gives results equiva- 
lent to those obtained by the previous method. Using the simpler 
method, we shall then study somewhat more complex interactions. 

We first present a brief outline of some results from the Appen- 
dix of Mathematical Biophysics (Rashevsky, 1938), using slightly 
different notation. Consider an axially symmetric cell of length ap- 
proximately 27, and of diameter approximately 27,. Let h; and D; 
be the permeability and diffusion coefficients of the substance 7. Simi- 
larly for the substance 7, let ¢); be the external concentration, ¢,;i the 
internal concentration at the ends, c;; the internal concentration at 


the circumference, and c¢; the average internal concentration. We 
have, then, for the net outflow over the ends (Rashevsky, 1938, p. 317) 


BB eat yaa et GANS art (ens leu) nye eee, (1) 


1 


and for the net radial outflow 
Aa, Di (¢; — evi) = ANT, 12 (Cri — Coi) i , (2) 


where we have taken 227%, and 427; Yo as the respective areas. 
If gq; is the average rate of consumption of substance 7 and 


= on 1, 772 is the approximate volume, the total outflow is —q@ V. 
Eliminating c.; and ¢; from equations (1) and (2), then adding the 
5 Sal BRAS ; 
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outflow from (1) and (2) and equating to —qiV , we have, after rear- 
ranging, 
Ci = Cop — Ai Gi (3) 


where 
2 (7, hi + Di) (%e hi + Di) 12 


~ 3h, Di 2710 i Di) + f2(M%2 i + Di) 


In the case of a sphere, 7; = 72 = 7, and we have 


(4) 


Aj 


Pl ae 9 F 
= (sara (5) 


BS 9 h; 


In section I of the preceding paper (Landahl, 1937), it was as- 
sumed that the glucose broke down to lactic acid at a rate propor- 
tional to its own concentration; that the lactic acid was resynthesized 
at a rate proportional to the square of its concentration; and that the 
lactic acid was oxidized at a rate proportional to the product of its 
concentration with that of oxygen. Instead of the differential equa- 
tions previously obtained, we have, upon using Rashevsky’s recent 
method, the algebraic equations 


Gp tH Cf} cys (6) 
G2 = Nf Cot, +- bc —ac,, (7) 
ian ROM (8) 
Substituting equation (3) into the above equations, we obtain 
Qi = @(Cor — Ai G1) — B (Coz — Ac G2)?, (9) 
G2 =NG,—h, (10) 
Gs = f (Co — Ae Qe) (Cos — As Qs). (11) 
Equations (9), (10) and (11) may be written 
Aq@.—Bqat+Cq—D=0, (12) 
da+@—nq=0, (18) 
Wes Ha4;—Kq:—Mqa+S=0, (14) 
A=bA*, H=f A. Az, 
B=20D 6 Ac, K=1-+ fee As, 
C= aA, 11, M = f es Ae, (15) 


D = 4 Co, — b €%p , S = f Coe Cns 
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Equations (12), (13) and (14) are the same as equations (22), 
(23) and (24) of the preceding paper (Landahl, 1937) so that the 
subsequent relations are formally alike. Equations (15) are the same 
as equations (25) of the previous paper if ¢, and e, are neglected, and 
if we use the above equation (4) for the definition of A;. In the case 
for which b = 0, the solution of equations (12), (18) and (14) for q, 
is given by equation (66) of the previous paper, in the form (¢ 3 = 2) 


Coo = BY (16) 
where y = q;/@s:, and q3: is the limiting value of q;, and 
1 180 
nce eee alle aia (17) 
A 
B = As qs1 = : [ (a Ay — 1) Coz - a Cor Ao] ° (18) 


N Ao (@ A, + 1) 
In equation (4) let us substitute 

m1=é2, Eslé (19) 
— is taken as greater than or equal to one since this is the more gen- 
eral occurrence. The results which follow will hold also for € < 1, 
but if not too small a fraction, i.e., for not too flat cells. We then ob- 
tain, omitting the index 7, 
— 2 (Er,h+D) (rr h+ DET 
| 3Dh 2E(Er,h + D)-+(mh-+D) 
Differentiating A with respect to €, we obtain 


ONS VEY 4h CR ot (2€r,h-+ D) (r2h-+D) 
fae DipiGEeDabDT 


Now for & = 1, we have A = A, of equation (5). For € = o, or 
Y, > > 712, we have 


A (20) 


Se ale 3 
eS Voor ie eae we 
Anno 2517 +5 Set ( ) 


where for 7, we use 7, the cylindrical radius. Since A increases with 
€ from A, to Bi so that the order of magnitude of A does not change 


with é , we may use the simpler expression A;, of equation (5) instead 
of A; of equation (4). Comparing the value of A (A;) of equation (5) 
with the definition of A in equation (25) of the previous method (Lan- 
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dahl, 1987) we see that inequalities corresponding to the previously 
derived equations (71) and (72) would lead to slightly larger values 
for h and values of D which are about three times as large. 

The value of A is determined both by the permeability, h, and the 
internal diffusion coefficient, D. The values of h and D are not deter- 
mined even approximately, without additional information. For this 
reasou it is often convenient to use A without reference to the per- 
meability and diffusion coefficients. A; measures the resistance to the 
flow of the substance 7 due to both the internal diffusion and the dif- 
fusion through the membrane. We may then refer to A as the total 
diffusion resistance. A has the dimension of time. The larger the A, 
the longer is the time required for a given.amount of substance to en- 
ter or leave the cell in a nonstationary diffusion state. 


II 


Let us now make an estimate of the value of the total diffusion 
resistance A, for lactic acid and the coefficient f [equation (8)] on 
the basis of data availabie for unfertilized Arbacia eggs. For low 
values of the external oxygen concentration c,;, the rate of oxygen 
consumption, q;, is limited by the value of c; and is proportional to 
it, while the value of the concentration of lactic acid ¢. is relatively 
constant. For large values of ¢ 3, the concentration of oxygen in the 
cell, ¢;, is high and stable, but the oxygen consumption is limited by 
the rate of production of lactic acid. (The latter holds if the glycoly- 
tic coefficient is very low; it is probably low enough in the case of un- 
fertilized Arbacia eggs. See later discussion.) Then near the point 
of maximum curvature of the y() curve, x = ¢ 3, both factors are 
of the same order of importance. If, at this point, Cc. = c*,, ¢; = C’s, 


and qs; = q*; approximately, then, as in equation (8), we may write 
d's.= fc, c's. (23) 


From the preceding discussion it follows that c*, and c*; are of the 
order of but less than the maximum values. Hence very roughly we 


may put ¢*, = } C2 max and c*; = 4 c*);. From equation (17) 


Taking C. maxx = 2 < 10+. gm./cc. as the approximate concentration of 
lactic acid in the unfertilized Arbacia eggs (Perlzweig and Barron, 
1928); and a = 12 * 10° gm./ce., cos = 6 < 10-7 gm./ec. and 
q*s = .85 X 10° gm./cc./min. = 6 & 10-8 gm. /cec./sec. as the approxi- 
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mate values obtained previously (Landahl, 1937) from the data by 
Tang (1931), we have 


f © 10° ce./gm./min. = 2 < 10° ec./gm./sec. 
(25) 
Ae © 8X 10? min. = 5 X 10? sec. 


In order that A, have the value estimated, we have, using the in- 
equalities of equations (71) and (72) (Landahl, 1937) 


Deo em min, = 250 10 em.2/sec. 
(26) 
hy © 2X 10° cm./min. = 3 & 1077 cm./sec. 


The concentration of lactic acid in the fertilized eggs is but slightly 
more than that in unfertilized eggs (Perlzweig and Barron, 1928). 
The oxygen consumption q;; is, however, four or five times as great. 
Since the other values are not very different, A, for fertilized Arbacia 
eggs is several times smaller than that for unfertilized eggs. 

Comparing the above values with that previously obtained for 
D; and hz , we see that they differ by a factor of 10°. On the basis of 
physical considerations one would not expect the diffusion coefficients 
for lactic acid and oxygen to differ by such an order of magnitude 
even though the former should drag larger molecules with it. On the 
basis of Stoke’s Law and the previous limit estimated, 


D, > 10-* em.?/min., a value D, > 10, 10° cm.?/min., 


would seem plausible. Hence the symbols are as chosen above in (26). 
The rough estimate gives a value for the permeability of lactic acid 
(lactate) of the same order of magnitude as that for glycerol ob- 
tained by other methods (Jacobs and Stewart, ’32). 

Equation (18), rewritten as 


1 Coz @ Cor 


yea eres el ae 


Op == 
gives the qz:(¢o2) relation and ” A, is the reciprocal of the initial 
slope. It would be of interest to compare the more direct evaluation 
of A. from the experimental relation q2:(¢o2) with the estimate made 
above. In the event that the experimental procedure to determine this 
relation is straightforward, it would furnish a method for the study 
of various factors on the membrane’s permeability to lactic acid. 
Among the factors of particular interest would be fertilization. How- 
ever, for cases of large glycolysis certain difficulties arise, to be dis- 
cussed later [equation (84), seq.]. 
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III 
Let us now consider a somewhat more general case in which the 
resynthesis occurs by two routes, one direct and proportional to b C5 
and the other coupled to the energy from the oxidation of lactic acid 
so that this part of the resynthesis is proportional to the oxygen con- 


sumed, v f G c;. Instead of equations (6), (7) and (8) we have 


% = 4C,—be»—vf ee, (28) 
Gp = Taf Cy Cp — @ Ca UC Vy Cree (29) 
ds = f Cs 6. (30) 


Substituting equation (3) into the above set of equations and re- 
arranging, we have, using equations (15), 


yHqq;+A G5 — (B+ »yM)q.— »(K — 1) 4 


(31) 
+Cqa—D+7S=0, 
a+ d@—nqa=0, (82) 
H @4;— KqQ:—-Mqa+S=0. (33) 


Multiplying equation (33) through by » and subtracting from (31), 
the latter may be written 


Aq@.—Ba+Cqa—D+7q=0, (34) 
Substituting for gi from equation (32) into (34) we obtain 
A q@,-— (B+ C)q.-+ (nC+7)qa,—D=0. (35) 


Comparing equation (35) with the equation resulting from the elimi- 
nation of q, from equations (12) and (13) of the preceding case, it 


is evident that n is replaced by (n -- a). Since (nC + yr) is constant 


and independent of c,; and since equations (14) and (33) are identi- 
cal, the form of the q;(¢ 3) relation, in particular, is unchanged. If 
the resynthesis is considered only coupled with oxidation so that b is 
taken equal to zero, we have as in equation (16) 


Cos = B Yy aa — ; (36) 
where now 
wy (a A, +1) 
y = ee eee 
f As[n(@ Ay+ 1) + 9] 
— As [(@ Ar + 1) Coo + & 9; Ag] 


Beas coc aot (Ce (38) 


(37) 
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IV 


The mechanisms of respiration thus far discussed are extremely 
simplified when one considers the actual processes which probably 
occur.* Let us now consider the effect of an intermediate product 
which does not leave the cell and which contains a radical which in 
no form leaves the cell. Consider the following reaction (Cameron, 
1930): 


a g 
glucose + phosphate wi hexosediphosphate > lactic acid + phosphate 
(1) (6) (7) Phen (2) (6) 
oxygen 


1|@) 


carbon dioxide + water 
(4) (5) 


Let us assume that the phosphate in no form leaves the cell, so that 
the total amount of phosphate per unit volume is a constant cp. For 
simplicity we shall not consider the reactions as reversible. Also, the 
equations for carbon dioxide and water shall be omitted. Then if the 
subscript 6 refers to the phosphate and 7 to the hexosediphosphate, 
and if m = 2 M,/M, and again n = 2 M./6 M,, where M; is the mole- 
cular weight of substance i, we have the following equations: 


= 46,6, (39) 
do = nf GC; — Cr, (40) 
ds = f C2 6s, (41) 
de =—; =406,6,—9¢, =0, (42) 
Cae mers (43) 


Consider, first, the anaerobic case so that C; = d3 = 0. Then solv- 
ing for ¢; in equation (43) and substituting into (42) and rearrang- 
ing, we have 


Caer oto (44). 


Substituting for ¢, the above value and using the relation ¢, = Co1 
— A, g;, we have, on clearing, 


am by ds — [emer tag erat g]ui+49 Cree =0. (45) 
* For summary and bibliography see Burk (1987). 
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For sufficiently small values of ¢o: , and hence qi, 


Ls a Cp A 
at acpA,+1 


o1 
and for sufficiently large values of Co1 , 1 approaches the limit 


en ai 
11 ™ 


If z = 4:/du , then equation (45) may be written 
n & 


where 


(46) 


(47) 


(48) 


(49) 


(50) 


Equation (48) is of the same form as equation (16) and may be 
illustrated by the curve in Figure 1, which is drawn to approximate 
the experimental values for fructose obtained graphically from data 


100 


Zz 
3 


_ LACTIC PRODUCTION RATIO 


PERCENT, GLUCOSE 


5 
PERCENT FRUCTOSE 


FIGURE 1 


The lactic acid production ratio, z = q,/q,,, is plotted against the concen- 
tration of fructose [°] (Boyland and Boyland, ’38) and glucose [.] (Warburg, 
26). For fructose, the value of q,, was taken as 6.9 X 10-6 gm./cc./sec. and for 
glucose 8.8 X 10-6 gm./cc./sec, Note that the abscissa scale is different for each. 
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by Boyland and Boyland, (1938). The experimental points for glu- 
cose are by Warburg et al. (Warburg, 1926, p. 121). Since the values 
of the lactic acid production are not given for higher sugar concen- 
trations, the highest value in each case was taken about ninety per 
cent of a theoretical maximum. The origin for the ordinates was 
taken, in both cases, as the value of the lactic acid production for zero 
concentration of hexose in the experiments by Boyland and Boyland. 
Since the number of experimental points in each case is so few and 
since the origin for the curve in one case is not certain, and the limits 
in neither case can be estimated with any accuracy, no attempt is made 
to evaluate any of the parameters. 

Because there are so few points, it is not possible, to say that the 
curves for glucose and fructose are significantly the same or different. 
But if the curves did coincide in the figure, then ys, = 20 yj and 
Cy, = 20,1. , where the subscripts refer to fructose and glucose, and the 
factor twenty arises because there happens to be a factor of twenty be- 
tween the two scales. Since the limiting values of the lactic acid prod- 
uct is nearly the same in each (Boyland and Boyland, 1938), then by 
equation (50) we have Aijr = 20 Ai or the total diffusion resistance 
for fructose is twenty times that for glucose. If we assume also that 
the secondary stage is to a large extent the same in both cases, so that 
cp is nearly the same in each, then, since m is the same for both, g 
must be the same for both to satisfy equation (50). In other words 
the second stage reaction goes at about the same rate in both. In this 
case one might then tentatively assume the secondary reactions to be 
largely equivalent. From equation (49) we have ay = 20 a;, , so that 
the rate of formation of the second product, hexosediphosphate, from 
glucose is twenty times the rate of its formation from fructose. How- 
ever, it could very well be that, at least in some instances, the two 
curves could not possibly be made to coincide by a linear change in 
the abscissa scale. In this case it would not be possible to compare 
the total diffusion resistance or other constants directly ; but it would 
be necessary to fit the curves empirically to evaluate 4 and ¢. Ai 
would then be given, but in general, an additional empirical relation 
would be required to evaluate a, g, and ¢p. 

In experiments on tissue extracts (Haarmann and Stratmann, 
1932) the results indicate an initial increase of glycolysis with sugar 
concentration and, in most cases, a decrease at higher values. The 
values are, however, quite erratic. In experiments on yeast fermenta- 
tion (Hopkins and Roberts, 1935) excellent curves are obtained. The 
initial rise is of the form of equation (48), but for rather high values 
of sugar concentration the fermentation rate drops very slowly. (The 
ratio of the initial slope to the final slope is approximately minus 
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fifty for glucose and minus one hundred for fructose.) The condi- 
tions in the two above cases are not directly applicable but are men- 
tioned because of their similarity to the case under discussion. 

If the supply of glucose is more than a few tenths per cent, 
dQ: == Gu. The aerobic case is then simplified and the result becomes 
a special case of section II of the preceding paper (Landahl, 1937). 
By analogy, we may write 


(51) 
where 
=_ (52) 
CRETE Fake 
B= As V1 = ~ (Cos aaa ay (53) 


V 

The glycolytic coefficient (Rashevsky, 1938, p. 105) §’ or Gly- 
kolyse/Atmung (Warburg, 1926, p. 126) is defined by the ratio of the 
number of lactic acid molecules produced by glycolysis which are not 
oxidized to the number of oxygen molecules consumed. Or we may 
say that 3 f’ is the ratio of the number of molecules of lactic acid gly- 
colyzed to the number oxidized. The symbol f’ is the same as the § 
used by Rashevsky (1938) and is not to be confused with the f used 
here. Then 


Boca Te PASI 0 2) ee (54) 


Solving for q. in equation (14), (33) or (41) and substituting 
above, 
jee Qs — f Coz (Cos — As Qs) 
32 f As Qs (Cos — Az Qs) 


Since y = q;/@s: and 6 = Az Qs: , we have, using equation (16), (36) 
or (51), 


(55) 


Cos —— Aa Das Cog Bells (56) 


We may then write for equation (55), using equations (87) and (38) 
quite arbitrarily, 


f= m(@A,+1) +» Mende (1—y) 


BAG te - Oem hI) ym 7) 
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from which it is evident that f’ decreases with ¢ 2, the external lac- 
tate concentration. In a great deal of experimental work (Warburg, 
1926) the value of ¢ 2 is negligible and for this condition 
y 1—y 
n(@ Ay + D! y } £58) 
For small values of the external oxygen concentration, ¢);, and 
hence small y , equation (58) may be written 


ek 
g==04 


gle SM: i GCytAs 1 
eae gee (aAgeeel) Cos bow) 
and for values of y not too far from unity. 
ave 1 
Lee Ter 7 Ure 20) 


For sufficiently large values of the external oxygen concentration Cp; , 
B’ approaches zero (or becomes negative if Co. 4 0). The coefficient 
B’ is inversely proportional to c,;. For a constant, not too high value 
of ¢;, 6’ increases with decreasing total diffusion resistance A, ; and 
to a lesser extent f’ increases as the external glucose concentration, 
Co,, increases; as the total diffusion resistance to oxygen, A;, in- 
creases, and as the total diffusion resistance to glucose, A, , decreases. 
Also f’ decreases with f and to a lesser extent increases with a . These 
results for the variation of f’ are essentially dependent only on the 
assumptions involved in obtaining equation (8), and hence apply to 
all the preceding cases. However, in the following case the expres- 
sion for the glycolytic coefficient, though similar in some respects, 
will differ in the respect that for large oxygen concentrations f’ may 
remain large and positive. 


VI 


Let us now consider a still more complex system and discuss 
some of its properties. Consider a system in which a sugar is glyco- 
lized through the hexosediphosphate stage into lactic acid. The mech- 
anism through which lactic acid is oxidized is complex. That in some 
instances intermediate stages are pyruvic acid, acetaldehyde, acetic 
acid and alcohol has been suggested (Cameron, 1930). Let us assume 
that the intermediate products do not appreciably leave the cell and 
consider a net reaction in which lactic acid must be in contact with 
a catalytic substance K, i.e., form a “complex,” in order that oxida- 
tion may occur.* We then have the following system: 

* It is also assumed that the catalytic substance does not leave the cell. How- 


ever, certain enzymes leave the cells in sufficient quantity so that the filtrate 
may be studied with respect to its metabolism (Rubenstein, 1933). 
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> hosphat ; lacti id + phosphate 
lucose + phosphate ~ hexosediphosphate ~ lactic acl p 
et) ‘ (6) i. (7) Se gee) (6) 
k K 
ae) 
(lactic acid) - (K) + oxygen ~ carbon dioxide + water + K 
(9) (3) (4) (5) 


The number under each substance is the subscript number, and cx is 
the total catalytic substance per unit volume. Considering only the 
forward reactions, we may write the following equations: 


Gge= Gicnewe (61) 
de=—7C, - NGC. (62) 
(ia CnC (63) 
Oe = Gr = UG Ce 0 Oe (64) 
(e=—O = a (65) 
Op = Cee cis (66) 
Cx =Ce+1ey. (67) 


As in section IV let us consider that the supply of glucose is am- 
ple so that the above set of equations reduces to 


G1 = Qi, (68) 
a, ——Ou +1 G3, (69) 
@=—ft;&=hG¢,, (70) 


together with equations (67) and (8). Eliminating ¢,, G , q@: and @> 
from equations (67), (68), (69) and (70) and using equation (3) 
we obtain after rearranging, 


[wk As+ fAs(1+ 2k Cx As) ] 92s — [ (Coo + Ae G11) 
xX (kL-+ kf ex As) + (Nk ex As +1)f Cos] Qs 
+k f €x (Coz + As Git) Cog = 0, (71) 
which can be given the same form as equation C10) ac HOr Co. oe 
43 = Qs1, and 


> Ke Cx (Co2 + As iz) 
Uk) == witele CL NEL Ae & mk PS 
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Equation (71) may then be written 


ol apes (73) 
where in this case 
eo ke U(Co2 + As Git) 
f (nk ex Ap 1)? oly 
(nkis,.+faA;tnkfe Ae Az) (Cog As Qu)ke 
= K ) -- qu) oS (75) 


f(wk ex Ag+ 1)? 


If there is an ample supply of the catalytic substance K present, 
we may let cx approach infinity in equation (71) and obtain 


[NM As As] Q?3 — [Az (Co2 +Av Gir) + 1 As Cos] Qs 

+ (Coz + Ag it) Cos = 0 (76) 

which may be written 
(1— y) (G3 —By) =0, (77) 
with £ defined by equation (53). An equation of the form of (76) is 
obtained if equation (70) is replaced by q; = f ¢ ¢; [ equation (8)], 
and if A, is very large, i,e., a is negligible. The latter restriction arises 
because of the assumption that the lactic acid catalyst does not leave 
the cell. In equation (77) the oxygen consumption increases linearly 
with oxygen concentration until it reaches its limit, after which it is 
independent of c,;. The equation (8) then corresponds to a special 


case of a reaction requiring a catalyst but in which the catalyst oc- 
curs in sufficient quantity not to be a limiting factor. 


VII 


The expression for the glycolytic coefficient for the case just pre- 
ceding, in which the sugar supply is large, is given by substituting 
the value of g. given by equation (69) into (54). We then obtain 


g pe? | Tih Sek (78) 

N Qs nN Qs N G31 Y 
For small values of the oxygen consumption, f’ is evidently large and 
unstable. Hence let us only investigate the variation of #’ for fairly 
large values of cs, where approximately ¢; = « o/ (A == etThen 


ae Qi —nk Cr Coz du b (79) 
ee Wh Guat Gita)? NT Oe ile Cx de 1) Os 
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and f’ decreases to a certain extent with the external concentrations 
of oxygen and decreases with the external lactate to a negative limit. 
Very frequently in experiments in vitro the external lactate is neg- 
ligible and the oxygen tension is sufficiently high so that the consump- 
tion of oxygen is independent of the latter. It is under these conditions 
that the glycolytic coefficient is generally obtained. For this case we 
have for the glycolytic coefficient 
; 1 

Mine eae Te (89) 
Now if for some reason certain cells undergo an irreversible change 
in which the available substance K is decreased, the total diffusion re- 
sistance to lactic acid A, is decreased, or the reaction rate k is de- 
creased, the coefficient f’ would increase in inverse proportion. 

Now for large ¢ 3, we readily find that 


q _ Wu — NK Cg Coe 


2 (Cos 223) ee ee * (3) 


For normal cells 3 6’ is much less than unity, seldom over one-tenth, 
so that n k cx As > > 1, and again neglecting co. , we have, with close 
approximation 


Ge aa: B Qit > (82) 


and only a small part of the sugar consumed is produced as lactate. 
For cancerous cells 3 f’ is of the order of one and often several times 
larger. For the more extreme cases we may write roughly 


; LE eens =— dp _ 


that is, nearly all the sugar consumed per unit time, qu, is trans- 
formed into lactic acid which leaves the cell. 

For normal cells we have for the limit of the oxygen consump- 
tion, from equations (72) and (50), 


GE = & Kegs ere! . (84) 


Equation (84) is the counterpart of equation (27) when ¢, is above 
the value at which it no longer appreciably affects the sugar consump- 
tion. The qs: (¢o2) relation in this case also gives A. as in equation 
(27 ). However, as we pass to cells with larger values of f’ the error 
in A, obtained from the qs: (¢2) relation of equation (84) or (27) be- 
comes larger, and for the more extreme cases 


Ast == K Cx (Coe + Ae Git) « (85) 
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The slope of the qs: (¢o2) tends to become independent of A, , and, in- 
stead of measuring the reciprocal of A, , it measures the product of 
the reaction rate k and the total catalytic substance ¢x . 

Equation (45) or (48) gives the relation between the aerobic 
production of lactic acid and external concentration of glucose. It 
would be of decided interest to compare accurately obtained experi- 
mental relations for both normal and cancerous cells of similar tis- 
sues under identical experimental conditions and over a wide range 
of values. It might then be possible to determine whether the differ- 
ence is primarily due to a change in permeability, or to the change of 
the chemical reaction rates a or g, or, finally to the catalytic action of 
cp. Additional knowledge would be obtained by the use of well 
known inhibitory agents. According to equations (49) and (50) the 
limit of the rate of glucose consumption, qi; , does not depend on the 
total diffusion resistance to glucose, A, , or on the reaction rate a, but 
the shape of the curve does. 

Similar parallel studies on the same kind of normal and patho- 
logical tissue for the determination of the relation between the aerobic 
production of lactic acid and the external lactate concentration and 
the relation between the oxygen consumption and oxygen concentra- 
tion, g3(¢o;), should lend some insight into the problem of malignancy. 
There are in the literature reports of a large number of experiments 
for some of the relations needed [for q;(¢ 3) relations see summary 
by Tang (1933)] and less on others. But the material used and the 
experimental conditions vary so much that any systematic study of 
the results from the point of view outlined above becomes difficult 
and of questionable value. It might be emphasized that, though the 
particular reaction mechanisms considered have some empirical justi- 
fication, they should be considered as approximations to certain ex- 
perimental conditions and are given to illustrate certain types of re- 
actions. It is of interest to note that the relation between oxygen con- 
sumption and oxygen concentration, q;(¢o3), is the same in all the 
cases discussed. [Cf. equations (40), (58) and (66) of the first paper 
and the above equations (16), (36), (51), (73). The previously ob- 
tained equations (40), (58) were not of the same analytic form but 
could not be distinguished graphically except perhaps under very 
favorable conditions.] If § in the above mentioned equations is fairly 
small we may write with very close approximation 


Co3 == Ae Sh (86) 
a) 
or 
31 Cos (87) 


Sem Apia gs 
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The equation in the form above has been empirically tested by Tang 
for a large number of cases (Tang, 1933). 


Summary 


The recent approximate method developed by N. Rashevsky is 
shown to lead to results equivalent to those obtained in the preceding 
paper on cell respiration. It is shown that, over a large range, many 
of the results are not dependent upon the cell shape. The results are 
simplified by the definition of A as the total diffusion resistance, the 
resistance to flow due to both the diffusion through the cell membrane 
and through the cell itself. Using information from various sources, 
an estimate is made of the total diffusion resistance to lactic acid, Az, 
and the reaction rate constant f for unfertilized Arbacia eggs. 

Several more general reactions are considered, such as the effect 
of the resynthesis being coupled to the energy of oxidation; also the 
effect of an intermediate stage in glycolysis, such as hexosediphos- 
phate; and the effect of a catalytic material required for the oxidation 
process. Expressions for the glycolytic coefficient are derived in terms 
of the parameters of the cell. An outline for a series of experiments 
is given, which, it is suggested, would be of decided value to various 
problems, especially to the problem of malignancy. 
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5 The approximation method, recently developed by N. Rashevsky, 
is applied to a problem of periodic cell reactions, studied exactly in 
a previous publication. The fundamental features of that previous 
investigation are thus obtained in a rather simple way. 


In a previous publication (Weinberg, 1938)*, we have demon- 
strated that in an idealized cell in which two coupled monomolecular 
reactions are proceeding, periodic variations in the concentrations of 
the reacting substances may be expected, at least if certain relations 
exist between the reaction constants and the diffusion coefficients of 
the substances involved. The method used was quite complicated, de- 
pending, as it did, upon the solution of a pair of partial differential 
equations; moreover the results were valid only for the special case in 
which s = ue where h; and D; are the permeability and diffusion co- 

a) v2 
efficient, respectively, of the ith reacting substance. 

Since the appearance of this article, Rashevsky (1938a) has de- 
veloped an approximate method for handling diffusion problems 
which simplifies enormously the mathematical detail required for the 
solution of these problems, and which enables one to treat cases which 
are much too complicated to be handled rigorously. In this note we 
shall show how this approximate scheme leads very directly to our 
previous results, and how it enables us to obtain, approximately, the 
periods in cases in which Ia ee 
1 4 

Let us again consider a spherical cell of radius 7), in which two 
substances whose average concentrations are ¢, and ¢, are being pro- 


duced according to the law 


Rate of production of ¢ = Au ¢:-+ Ar C2 
(1) 


Rate of production of 6: = As: 6, + As. 
Ag will be seen presently, the introduction of the average concentra- 


* Referred to as (A). 
“21925 
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tions reduces the partial differential equations previously obtained 
to ordinary ones. a 
If c,* is the concentration of the first substance just inside the cell 
membrane, it is clear that the average concentration gradient of this 
Es Be 


ee Css 
substance is ————— 
To 


of material passing through a square centimeter of surface at the cell 
wall, 


. Fick’s law then gives for the average amount 


D,(¢, — ¢1") 
To 


2 


and since the total area of the cell surface is 42 7, , the total amount 
passing through the cell wall is 42 r) D, (¢, — ¢,*). Denoting the aver- 
age rate of production of ¢, per cubic centimeter by ¢, (so that the 


total amount produced per second throughout the cell is 3 nN 4 a); 


we find from a simple “material balance’’ consideration that 


4 es = 4 dc, 
tod — 42% D, (¢: — 61") fe is ao 


(2) 


A similar equation holds for ec, with subscripts 1 replaced by 2 every- 
where. 


To determine c,*, we have the boundary condition 
DMG sos 
hy (¢1* — ¢o1) = steal (3) 


Yo 


where ¢o, is the external concentration of ¢, and is assumed to be con- 
stant. Solving for c,*, 


Die, +% In Con 


PREMECSUUES oe 


If we substitute in equation (2) the values for q; and c,* given in 


es (1) and (4), we finally arrive at the pair of ordinary equa- 
ons 


dt, 3h,D ra tt 3h, D 

Stn Ah noe le Dee ie. 

hme Ae MERE MR abe korn area 

dc: 3 Ny D. ” 
2 Fe ere | 

<a =e Ans Cy A ae eK LE Se 

Eo 0 et Aw Fr hac) 2 eae eee 
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The solution of a set of equations of this sort is perfectly straight- 


forward (see, e.g., Rashevsky, 1938b). The concentrations ¢, and ¢c 
are found to be essentially exponential functions of the time of the 
form 


Ci ay eo. = Bi 
where a; and /; are constants and A is given by 


m?, D,; +m’, D, 


i 4 (Ay + Aza) a D 


+ V4 (Au — Ane) — (m*, Di — mz D2) ]? + Ay An. (6) 


3h; 


1) (Di + 1 hi) 
Equation (6), which defines the period of oscillation of the con- 
centrations, is entirely general. If the simplifying assumption intro- 


The symbol m?; stands for the quantity 


h 
duced in (A) is made again, i.e.. — = Ia 


Dz Al we find that m?, = m?, 


= m?, and so 


A= 4(An+ den) —2 (D+ Ds) + (Di —Ds) 


xq) EE An — An) /2(D, — DPF Ae An/ (D1 —D2)* (D 


which is formally the same as equation (34) (A). 
In the rigorous treatment, the quantities m? were determined as 
solutions of the equation 


EY OU el eerren wen (8) 
whereas in this case 
3 hy 
p ariehchtens DAW io (9) 
in 1) (D, + 7, h;) 


That the first non-zero root of (8) (corresponding to the first har- 
monic) is of the same order as the m? given in (9) may be seen by 
considering two limiting cases. First, when D, > > 7o h,, as seen in 
(A), the first non-zero root of (8) is 


3h, 
D, 1 


mM = 


? 
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which result also follows directly from (9); and second, when 
D, < < 1% h,, (8) becomes approximately 


tan m7’, = 0 


and so 
n? 
Mm=—, wWwoe—, 
To To 
while (9) gives in this case 
3 
Mm? = ar 
ro 


which is about 1/3 the above value. For intermediate values of D we 
should expect better agreement between the two results, since the two 
are identical for large D. 

We have thus shown that the period of oscillation of concentra- 
tions predicted by the use of this simplified method agrees very well 
with the period of the first harmonic of the rigorous solution (A), 
and thus checks the correctness of the previous solution. Furthermore, 
the approximate method has enabled us (Eq. 6) to determine the 
period of oscillation when is we es , a problem which could not be 
solved rigorously. 

The writer is glad to acknowledge his indebtedness to Professor 
N. Rashevsky for having pointed out to him the applicability of this 
method to the case at hand. 

The preceding results may be improved by doubling the value 
used here for the average concentration gradient, in accordance with 
Rashevsky’s considerations in a more recent paper (Bull. Math. Bio- 
physics, 1, 23-30, 1939). This modification changes the definition of 
m? so that (9) becomes m? = 6h,/7,(2D, + 7 h,) and it makes the nu- 
merical agreement with the rigorous result better, although it leaves 
the general] result unchanged. 
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THE MECHANISM OF CELL DIVISION 
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The recently developed approximation method for treating prob- 
lems of cell biophysics is generalized and corrected in some points. 
The equations of elongation of a dividing cell of any shape are given 
for the most general case of finite permeability as well as of finite 
internal and external diffusion coefficients. 


In the Appendix to our book (Rashevsky, 1938a) we have devel- 
oped a general approximation method for treating problems of cell 
deformation under the influence of diffusion. The method while only 
approximate has the advantage of generality, applying to any elon- 
gated shapes without restrictions to an analytically definable boun- 
dary surface. At the same time it leads to a very simple and short 
way of treating such problems, for which an exact mathematical treat- 
ment requires a rather involved procedures. H. D. Landahl (1939) has 
applied this method to the problem of cell respiration and has ob- 
tained in a much simpler way the same results which he found pre- 
viously (Landahl, (1937) by the more elaborate old method. Also A. 
Weinberg (1939) obtains by this method essentially the same results 
concerning periodical reactions, as those which are given by the much 
more elaborate exact method (Weinberg, 1938). 

A comparison of the equations for the rate of elongation of a di- 
viding cell, obtained in the Appendix (Rashevsky 1938a), to those 
obtained in a particular case exactly has been made by G. Young 
(1939), who worked out the exact expression for an ellipsoid of revo- 
lution. This comparison again shows the applicability of the approxi- 
mate method. The differences found are of no practical consequence, 
since no dividing cell has a regular shape, and the application to such 
a cell of formulae which hold exactly for an ellipsoid only would be 
likely to involve as large an error as that by which the approximate 
equations differ from the exact ones. 

In the original exposition of the method given in the Appendix 
to our book (1938a) we have for sake of simplicity considered the 
elongation of the cell only in a very special case, namely when both 
the permeability h and the external diffusion coefficient D. are in- 
finite, although the diffusion problem was treated for the more general 
ease of finite h and D; (Rashevsky, 1938a, p. 317). As has been 
shown by G. Young (1939), in the case of h = D, = o however, we 


29 35 
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should obtain a zero elongation, since c is constant along the surface. 
The error committed by us lies in the discussion on p. 308 (Rashev- 
sky, 1938a), where the force appears as a small difference of two 
quantities. The use of Betti’s theorem (Love, 1906, p. L7Z; young, 
1939), combined with our approximate method for determining the 
distribution of concentrations, leads to a much more straightforward 
derivation of the elongation equation and avoids the above mentioned 
doubtful argument on p. 308 (Rashevsky, 1938a). 

Betti’s theorem, applied to a plastic body (Young, 1939), gives 
for the average rate of elongation of a body of any shape of volume V 
in the direction of the z-axis: 


ak OA es 


1 
oe oe nT s,v 555 [z Z— BUR Sarton kA Rares 7-8) [Zy z 


=F (Xe 0 Yon U lis 


where X, Y and Z are the components of the force per unit volume 
and X,, Y, and Z, are the components of any surface force. The 
triple integral is extended-over the whole volume of the body, the 
double integral over the surface. 

If X, Y, Z are derived from a potential, so that (Rashevsky, 
1938a, Chapt. VII, VIII, IX): 


___ Flu 
ae M 
then as is shown by G. Young (1939), the volume integral transforms 


into a surface integral of the form: 


RT u 
— ap SS (cos (», 2) — $ [a cos (vx) + y cos(y, y)]} dS. 


grad ¢ 


Therefore we have for the average rate of elongation for a body of 
any shape exactly: 


—— T (1) 


where, taking z as the direction of elongation: 


; J = Sf {zcos(vz) —4 [x cos(vx) + y cos(vy)]} dS (2) 
an 


T=. (Sil Zve—4 (Xen Vy) as © (3) 


X, Yy Zy being the components of the force kit (cj; — c.) at the 


membrane. For h = D,. = w, (¢. — ci) = 0 and ¢ is constant along 
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the surface, which makes both J and T vanish. It is therefore indi- 
cated to study the more general case, where D;, h and D, are all 
finite. We shall then substitute the approximate value of c at the sur- 
face into (1), and thus obtain an approximate expression for the rate 
of elongation. 

In case of a finite D, the average concentration at the membrane 
just outside the cell is no longer equal to c. We shall denote it by 
c’, for the “ends” of the oblong cell and by c’, for its “sides”. The 
average diffusion flows just outside the cell membrane per cm? are 
proportional to c’,; — c) and c’, — ¢ respectively, and to D.. Denot- 
ing the coefficient of proportionately by 1/6 we have for these flows: 

a (c's — Co) and (c's — Co)- (4) 

re) 
6 has the dimension of a length and represents a sort of average dis- 
tance from the cell, at which the external concentration practically 
becomes equal to ¢ . 6 is of the order of magnitude of the cell size, 
but the exact relationship will vary with the details of the cell forms. 

At this point we shall also make a slight improvement of the treat- 
ment originally presented (Rashevsky, 1938, Appendix). For a not 
too irregularly varying concentration inside of the cell, the former 
will have an extremum in the middle, while the average value c will 
be about half-way between the middle and the periphery. Therefore 
in expressions (48) of the Appendix of our book it is better to use 
4a and 4b instead of a and b. This does not alter the shape of the 
final equations, but improves the numerical values of the coefficients. 
In fact, if we calculate ¢ for a spherical cell by averaging the exact 
expression (Rashevsky, 1938a, Chapt. II, equ. 22), we shall find, for 
D,. = ow, an expression of the same form as equ. (54) of the Appen- 
dix, except that the coefficient of q a?/D; turns out to be 1/15 instead 
of 2/9 while the coefficient of g a/h is 2/9 instead of 1/3. The above 
proposed improvement makes this coefficient in the approximate for- 
mula equal to 1/9 leaving the coefficient of ag/h equal to 2/9, as 
against 1/3 in the exact expression. 

For the stationary state our fundamental equations become now, 
by the same line of argument as used before (Rashevsky, 1938a, Ap- 
pendix) 


D,(¢ — C1) = $ah(c,—Cc',) (5) 
D;(¢ — ¢2) = 4 bh(c, — €’2) (6) 
2D; - 


Zee Ta AT ae) (7) 
a f) 
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at (OA) — (C's — Co) 


D; es (Cm Fer) ea (Oc) 5 aba 


Equ. (5) and (7) give, eliminating c’:: 

Pe (DD. — dD ;h) ef ahD.Cy 
2D;D-. + 26D;h + ahD, ; 

Similarly from (6) and (8): 


i 2; (D;,D. a dD;h) C —- ohD Co 
(2 07D, 26 Deeb ne 


C—— 


(8) 


(9) 


(10) 


(11) 


For D. = o, (10) and (11) are of the same form as (52) of the Ap- 
pendix, except for the correction by the factor 1/2, as should be the 


case. 
From (10) and (11) we have: 


vo, SOT Dexia 28 Deane De 
- bhD.(€ — Co) 
C— = 


2D,D. + 26Dih + BAD, © 
Substituting (12) into (9) we find: 


7 abq 
ham Rite, Des 


(2D;:D. + 26D;h + ahD.) (2D;D.+ 26D;h + bhD.) 
2(2D;D. + 26D,h + ahD.)a-+ (2D;D.+ 26Dih-+ bhD.)b 
For a sphere (a = b) this gives: 


2 aq 1 a’q 2 adq 
cori tas Sia aes 
with 6 ~ a this is of the right form. 
From (10) and (11) we have: 
Co — Cy 
2(DiD. + 6Dih) hD, (a—b) (c—ey) 
~ (2D;D, + 28D;h + ahD,) (2D:Ds-- 25D OAD.) 


(12) 


(18) 


(14) 


(15) 
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Introducing (13) into (15): 


abq 
ee ag 


(16) 
2(D.-+ 6h) (a —b) 


2(2DiD. + 26Dih + ahD.) a+ (2D,D.-+ 26D;h + bhD.) b © 
For h = o, (15) and (16) reduce to 
26D:D- 


ee Se (On ane OaD aL DD as mac te (15’) 
and 
a abq 26 (a— b) J 
Pik talon Disa eI aa sD, oD) 6 nto? 
From (5) and (12) we have: 
Cc, — Cy Me SOD og (c—c ) 
2D;D. + 26Dih + ahD-. - 
(17) 
Ate 2D,D. = 
as ee OT Os NEADS Cay eens 


For h = «, c, — ¢; = c — c, = 0 as should be the case. 

Now we may calculate J and 7. On the “sides” cos (»y, z) = 0, 
while x cos (v, x) + y cos (vy, y) = b. Hence for the sides the con- 
tribution to J is 


—C<$bX 2ab K 2a = — Ven. (18) 
In a similar way we find that the contribution of the “ends” to J is 
Vc,; hence 
J = V(¢,—c.). (19) 
In calculating T we have 


RT 
Xe + Yvy = = Oka — ©) 


while 
RT 

LvZ = 71 a(¢, —_— Ca)s 
Therefore 

RT ‘ ’ 

eM hee 7 gel OU Nee Ae A ect 1 (20) 

Introducing (19) and (20) into (1) and substituting into that (15) 
and (17) we find: 
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1lda_ 2kT 
adt 3My 
(u—1)D,D2h + 6DiD hip (21) 
x | ! 
“S (2D;De + 26Dih + ahD-,) (2D; De +- 25D;h + bhD-) 


x (a —b) (e—). 


For a = b the force and da/dt vanish, as should be the case. 

For ¢ > 0,¢—& > 0; alsoa—b > 0, (u—1) < 0 always. 
Hence, for sufficiently small 6 ~ a ~ b, da/dt < 0 for g > 0. If how- 
ever 6 © a~ b is sufficienty large, da/dt >0 for g > 0. This agrees 
with Chapt. IX, MB. To insure positivity of da/dt we must have 


Picteealenr T0010 cm. (22) 
For h = o, D;, D- finite, (21) reduces to 
lda 2RTu 6D;D- 


adi 3My (23D, 0D) @s0;-b6D).) a gees 


For h = wo, D;, D. finite, (18) reduces to 


- abq (26D; + aD.) (26D; + bD-) (24) 


°° 3D,D. (26D; + @D.)a-+ (26D: + bD.)b 
Substituting (24) into (23) we have: 
lda 2RTabqu (a— b)6 
adi 3 My 25D, eD.)a+ GOD. oD)6 ° 
Putting 6 = b, (25) becomes: 
Lda _2RT (dra Dodi iy aul ra 


a dt 9Myn 2(2bD;+ aD.)a+ (2D; + D-)b? 
For a >> b this gives: 


lda RT qu : 
adt 9MnD,- 


a form to which equ. (22) of the Appendix (Rashevsky, 1938a) re- 
duces also fora >> b. 

_ G, Young, 1939, p. 31) makes a comparison of equ. (26) with his 
equ. (50), which holds exactly for an ellipsoid. The two equations 
are in substantial agreement, especially over the range of elongations 
which alone comes into consideration practically. 
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Since no way has been found yet to integrate exactly the diffu- 
sion equation for any other shape than the sphere for the case of a 
finite h , equ. (21) may be used as an approximation in this general 
case. 

Adding to the right-hand side of (25) the effect of the surface 
tension (Rashevsky, 1938a, p. 311), we obtain 


Lda  2RT ab qu (a—b)6 
a dt 9My  2(26D; + @D.)a-+ (26D; + bD.)b 
ea) 
mar Sta) ste 


Requiring da/dt = 0, and making a = b = 6, we find from (27) 
as the critical radius of a sphere below which for any elongation 
da/dt < 0, that is below which it will not elongate 


= Eee, + Dy 


— oRT iq ead 
This is in better agreement with equ. (34) Chapt. VIII (Rashevsky, 
1938) than equ. (34) of the Appendix, because in (34) Chapt. VIII 
D must become infinite with D. = o. Equ. (28) requires that the 
size of a cell should vary as q-”?. Since the “effective” q is propor- 
tional to the oxygen consumption (Rashevsky, 1938, Chapt. xX), we 
can check the above relation experimentally. A survey of available 
data shows that the required relation is approximately satisfied over 
a rather wide range of values of g. (Rashevsky, 1938c). 

As has been pointed out, for sufficiently large elongations equ. 
(26) reduces to that used before (Rashevsky, 1938a, Appendix). The 
latter was compared with experimental data on the elongation of a 
dividing nucleus, observed by K. Belar, and found to represent those 
data satisfactorily (Rashevsky, 1938c). : 

In the light of the present discussion the results of our recent 
paper (Rashevsky, 1938b) concerning cellular forms would have to 
be modified as to details, although the general methodological results 
remain the same. 

The author is indebted to Mr. H. D. Landahl for checking the 
algebraic calculations. 
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_ After deriving some auxiliary equations for the average elonga- 
tion of a viscous body under the action of forces derived from a po- 
tential, the diffusion problem for an ellipsoidal cell with a constant 
rate of reaction is solved for the case of an infinite permeability. 
The equation of elongation of such a cell under the influence of diffu- 
sion forces is derived, and compared with the approximate expres- 
sion obtained by N. Rashevsky for any kind of oblong cell. The two 
equations are in fair agreement. Effects of constant and variable 
surface tension are studied. 


Average rates of strain in a plastic body: 

One of the basic problems in the theory of elasticity is to deter- 
mine the deformation of an elastic body under the action of specified 
volume and surface forces. The general problem is difficult and in 
practice solutions are known only for some particular cases. How- 
ever, certain quantities connected with the deformation can be evalu- 
ated readily by means of results due to Betti. For example if w(xyz) 
denotes the displacement along the z axis of the point of the body 
originally at (xyz), then the average value throughout the body of 
6w/6éz is given by a simple integral involving the applied forces. The 
quantity 6w/6z is the relative elongation, or strain, in the z direction. 

Now in the case of incompressible material (Nadai, 1937, p. 421; 
1928, p. 467) (Goodier, 1936) there is an analogy between elastic and 
viscous bodies, in that strains in the first correspond to rates of strain 
in the second. If W (xyz) is the velocity of displacement in the z direc- 
tion, then the rate of strain in that direction is given by e = 6W/6z. 
The mean rate of strain in the z direction, defined as 


e=<Sffedv, (1) 
ean be found by employing Betti’s formula (Love, p. 172) in connec- 
tion with the above mentioned analogy. The result is 

aie (eZ—4(aX+yY¥)} de (3) 
Bnvlt sf (eZ —t(axX+yV dS. 
where V is the volume of the body, 7 is its coefficient of viscosity, 


ih: 4 Peet 


é.. = 
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X Y Z are the components of applied force per unit of volume, and 
X’ Y’' Z’ are the components of applied force per unit area on the body 
surface. 

Kinematically the quantity é is interpreted as follows: If e had 
the constant value é throughout the body, then the length of z com- 
ponent of a line joining any 2 points of the body would increase ac- 
cording to 


pence e One (3) 


This expression also holds if the body is a cylinder and e is a function 
of z alone. More generally, the expression 


ve=SSSeav (4) 


may be transformed (Phillips, p. 70) into 

ve=Sf Weos(nz)dS.- (5) 
Thus é is determined by the shape of the body and the z velocity of 
each part of its surface; no knowledge of the motion in the body in- 
terior is required. 


If the body forces derive from a potential the volume integral in 
(3) may be further simplied. Letting (X Y Z) = grad © the integral 


6® 
has terms J, = fff za , etc. Integrating by parts with respect 
to z gives 
Jz= Sf @zcos(nz)ds—fffodyv. (6) 
8 v 
In the same way 
Je= Sf g¢ucos(nx)ds—fffedrv, (7) 


and similarly for J,. Hence the total integral is 
J=J,—4(J2+J,) = Sf © [z cos (n 2) — 4{x cos(y x) 


+ycos (ny)}] dS. (8) 


This is seen, upon comparing with the surface integral term in 
(2), to be the same as if a pressure — ® were acting normally over 
the body surface. This is related to the fact (Love, p. 217) that the 
effect of a system of volume forces acting upon an elastic body can 
always be duplicated by a set of surface forces. . 
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If is constant on S then (8) becomes 


J=@ SJ [z cos (nz) —4{x cos(nx) + ycos(ny)}] dS- (9) 


which vanishes identically. This may be seen by considering that 
SS 2cos(nz) dS = V, and hence that the integral (9) is essentially 
V—i(V+ VD). 

Hence a constant surface pressure, or a set of volume forces de- 
riving from a potential which is constant on S, gives zero average 
elongation in any direction. 


Diffusion in an ellipsoid: 


Previous studies of diffusion processes in cells have for the most 
part been confined to spherical or nearly spherical shapes. While this 
may suffice for many cells in their normal or resting states, in which 
they are nearly spherical, it is not adequate for other important cases. 
In particular the phenomena of cell division involves changes in shape 
which carry the cell well away from spherical form. 

N. Rashevsky (1938, Appendix) has recently given an approxi- 
mate method of evaluating a number of quantities concerned in the 
diffusion process for a rather wide class of elongated cell shapes, such 
as cylinders and ellipsoids; and it is of interest to have some more 
exact results available for comparison. 

The next simplest shape after a sphere is an ellipsoid of revolu- 
tion, and this will be treated in what follows, the only approximation 
being in the assumption that the cell membrane offers no resistance to 
the passage of the diffusing substance. The results thus obtained are 
found to agree fairly well with those of Rashevsky and hence afford 
a certain degree of substantiation of his general method. 

Attention will be restricted to a prolate (elongated) ellipsoid of 
revolution inside of which the dissolved substance is being produced 
at the constant rate q per unit of volume. Let D represent the co- 
efficient of diffusion and c, the concentration of the dissolved sub- 
stance at great distances, and let the subscripts “¢” and “e” refer re- 
spectively to the regions internal and external to the ellipsoid. De- 
note the ellipsoid surface by S. Then the equations to be satisfied are 


V2¢,=2k, 
(10) 
2k = —q/D;, 
V?ce=0, 
(11) 
C.-C ato, 


i=, 00 (12) 
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Dey ees (13) 

Om) 61 
Taking the z direction along the axis of rotation a suitable sys- 
tem of co-ordinates (MacRobert, p. 198; Young and Reiner, 1937) is 


given by 


x= Vr? —f’ sin é@cos ¢, 


y= Vr —f? sinésing, (14) 
B= PCO - 


The foci of the ellipsoid are at (0, 0, + f), and the semi-major and 
semi-minor axes will be denoted respectively by a and b = \/a? — f?. 
A special solution of (10) is obviously 
ee — PCOS a, (15) 
to which may be added any solution of V? c; = 0. Thus 
rT 


pee 6). (16) 


¢,=kr? cos? 6+ 3 Ay Pn ( 
0 


Similarly the solution of (11) is 


Ce = Co-+ S Bu Qu(Z) Pp (C08 6) « (17) 
0 
Herein P, and Q, are the Legendre functions of the first and second 
kind. Since c¢, and c; are even functions of z it is seen at once that 
A, = B, = 0 for n odd. 
Expression (15) may be written as 


ke 795 P, (cos 6) +5 Pa (cos 6)]- (18) 


The boundary condition (12) at the ellipsoid surface (r = a) then 
gives, upon equating coefficients of P,(cos 6) , 


k 2 
<= + As =Co4 Bo Qo(S) » 


f 
2k a? a a 
A, P;(—) = 2 Yel) > 
3 - Ps (>) B. Q.(>) (19) 
Ay P.(=) hp Q.(F) forn > 2. 


The boundary condition (18) is equivalent to 
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Di =D, ate =a, (20) 
which gives 
shah = TB. Qo(2), 
shaft As Ps() =5°B,Q(4), (21) 
A, P’s(=) os a Be Qn) forn > 2. 


From the last relations of (19) and (21) it follows that in gen- 
eral A, = B, = 0 for n > 2. Solving for the remaining coefficients 
gives 


i Line? DOS 0 
Ay = %— ska Bean O° 
Deen DX tags 
Jassie aya 
A eet Oe ce 
* . Bie D Ps Qy—-DFaQs 
2 D; 
Bemigh tla -p-o,— DP.” 


where the omitted arguments of the P and Q functions are all a/ fe 
Thus the solution of the diffusion problem is 
r 
Ce = Co-+ Bo Qo G) + Bs Q:() Ps (08 6) ; 
(23) 


c; = kr? cos? 6-+ A,-+ Ap P, (5) P, (008 6) ; 


with the coefficients given by (22). 


Average concentration inside cell: 


One of the quantities considered by Rashevsky in his treatment 
of non-spherical cells is the average concentration of the diffusing 
substance inside the cell. For purposes of comparison this quantity 
will now be evaluated for an ellipsoid; namely 
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Bnel igt 24) 
Gee al (eC ( 
Ue 


In ellipsoidal co-ordinates the volume element is 
dv = (7?— f? cos? 6) sinodrdédd¢, (25) 
and the range of integration is on ¢ from 0 to 2 a, on @ from 0 toa, 
and on r from f toa. Omitting details of the integration, the result is 


@= Ay s (kat— 42). (26) 


For the case D, = o (22) reduces to 


1 
Ag Co ae 


24 a 
Ase ee (27) 


Bo = Be 08 
Substituting in (26) and simplifying gives 


am q a? b2 


CSS Da eee ey 


which should be compared with Rashevsky’s expression (Rashevsky, 
1938, equ. 24), in which we must make D, = ~, 
ws az b2 
C— 165 =f eee (29) 
3 i 202 
His a and 6 denote the half-length and the radius of a “best-fitting” 
cylinder of the same volume as the body. If these are chosen so that 
a:b=a: b then (29): becomes 
a? b2 


06 = 0.26 Foe (30) 


which is the same order of magnitude as (28) but about 25% larger. 
Mechanical effect of volume diffusion forces: 


A substance diffusing through a stationary material eerie on 
the latter a drag force F proportional to the concentration gradient, 
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fa Ve 
(31) 
=—RTu/M, 


where M is the molecular weight of the diffusing substance, and F is 
the Boltzman gas constant. To conform with the particular picture 
of the structures of the cell used by Rashevsky (1938; chap. viii), we 
introduce the indefinite factor « © 1. Further detailed studies of the 
effect of cytoplasmic structure on the mechanical properties of the 
cell may either justify the omission of u, or a general modification 
of (31). 

If then diffusion is taking place in a body we should expect the 
latter to be deformed by these forces. Regarding a cell as a viscous 
body (Rashevsky, 1938, p. 130) we can from the above results deter- 
mine the average rate of elongation of an ellipsoidal cell under the in- 
fluence of the internal diffusion forces. Since (2) is linear in the force 
terms, the contributions to é of any other kind of forces can be cal- 
culated separately and then added thereto. 

For a system symmetrical about the z axis the x and y parts of 
the integral in (2) are equal, so that (8) becomes 


J=Sffsf@Z—«ceX)dv= SS oz cos (nz) —x cos(nx)] dS .(82) 
From (23) it follows that on the ellipsoidal surface S 


c; = constant + A cos? 6, (33) 
where 
3 a 
A= kat 5A: Pal). (34) 


In ellipsoidal co-ordinates (Young and Reiner, 1937) 
cos(nz)dS = b? sinécosadgdd¢ 
cos(na)dS =ab sin? 6cos¢dédd. oo, 


Writing ¢ = a ¢;, and using (14), (33), and (35) in (32), we find 
upon evaluating the integral, in which the range of integration is from 
0 to220n¢ and from 0 tozoné, 


_ 8na 


is 15 


ab?A ; (36) 


and thus the average rate of elongation along the axis of revolution 
which would be produced by the volume diffusion forces is, according 


to (2), 
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2aA 
oo ies 387 
C= fo nw 157 (37) 
From (22) and (34) the quantity A is found to be 
Po 20 

fry ee f a ; (38) 

D; P's — De P2 ~ 

Qe 


where as before the P and Q functions have argument « = a/f. From 
the familiar expression 
P.(“) =4(32?—1), (39) 


the numerator of the fraction in (38) is simply f ; while the denomi- 
nator may be simplified by using the relation (MacRoberts, p. 114) 


P, Q'. — P’, Q2 = —1/(x? — 1). .(40) 
The result is 
‘Ar ___ sisi WES ate aa (41) 
3— (D; —D; OO 
j ‘alt aig Coe= 5 ING 
For a nearly spherical shape x becomes large and Q, becomes 
26h 
liar id al gaye 
Q(2) = Bote (42) 
and in this case A reduces to 
2 D; 
SSeS 2 SO 
= ae 2D,+ 3D. ee 


For elongated shapes x tends to unity, and Q:(a) becomes infinite. 
However it can be shown that 


Lim (2? —1)Q2(x) = 0; (44) 


so that the second term predominates in the denominator of (41). 
Then 


A=kafot (a? —1)Q,(2). (45) 
Explicitly (MacRobert, p. 114) 
Q(x) = H(8 2? —1) log ZTE _ 6 9). (46) 


Near x = 1 the first term of this predominates, so that effectively 
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32°—1 e+1i 38a?—1 (~-+ 1)? 
a ae eil tae ae 


Sees | a 
Q: dae as od Tie ML Psa 
(47) 
3 4v?—1 
sere [log (a + 1)? — log(«?— 1)]- 
Here the second term predominates leaving 
Q?(x) = — 4 log(a2* —1); (48) 
and thus we find 
ate, Bt ee b? 
ate Dc log =? (49) 
which vanishes as b goes to zero. 
From (10), (31), (87), and (41) it follows that 
: a Sa he 2. h? 
prementy Teer Pod (50) 
5s M y 2 D; + ™m D; 
where 
2 
(51) 


Ue apare(as Set) Op GTO 


By use of formula (46) this may be evaluated for different values of 
a/b. The result is given in the first two columns of the following 
table: 


a/b m te m 

22.4 142 25.4 277 

10.0 42 12.3 58 
Y fal 26 9.3 30 
5.1 16 7.15 16 
3.6 10 6.1 9 
2.4 6 4.3 4.7 
1.8 4.6 3.8 3 
1.34 3.9 3.5 Paral 
1.15 3.3 ai) 1.8 
1.07 3.06 3.23 1.65 
1.00 3.0 ee 1.57 


The last entry, for a = b, comes from (43). As the body elon- 
gates the effect of the diffusion forces tends to depend more and more 
upon the value of D, alone, the D,; term gradually dropping out. 

It may be noted that a fair approximation, up to a equal to about 
10 b, is given by taking m = 3 a/b. 

This should be compared with Rashevsky’s approximate expres- 


sion No. 26 which, witha: 6 =a: 0,is 
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RG a—b)a.b? 
tae — 0.17 RRs 6 ie 7 Be) 
adt My D,2b(2a+ 6) + D.(2 a + 6’) 
In effecting the comparison we write (52) in the form 
Dies Se 
Lda _2RPae oa o 
a dt 45 My tD,+m ID). 
and calculate ¢* and m’ for various values of a/b. Explicitly, 
t* = .5384 (a+ b) (2a-+ b)/adb, (54) 
and 
m’' = .262 (2a? -+ b?) (a+ b)/ab?. (55) 


Values of these quantities are given in the 3rd and 4th columns of the 
above table. If D, > > Dj; then expressions (50) and (53) are in- 
versely proportional to m and to m’ ; and Rashevsky’s value for é 
varies from 190% of (50) for nearly spherical shapes down to 70% 
fora = 10b. If D; > > D,. the expressions are inversely propor- 
tional to 2 and to t*, so that (53) is somewhat less than (50). If 
D,; = D, the agreement is still better for shapes not too far from 
spherical, as is seen by comparing m + 2 with m’-+ @t. 


Effect of constant surface tension: 


Another type of force involved is that of surface tension. A sur- 
face tension constant over the cell surface tends to keep it rounded 
up, and so opposes any elongation. We now calculate this effect for 
an ellipsoid of revolution. 

The principal radii of curvature of such a surface are 

(a? — f? cos? 6) 8/2 


R, = PS PP me sy (56) 


which is the curvature of a meridian section, and 


Re 3 (a2 — f? cos? 6) ¥. (57) 


, i 1 
The mean curvature is K = RP -+- Ro’ and at each point of the body 
1 2 


surface a surface tension y exerts an inward normal pressure equal to 
y K. Thus the surface stress in the z direction is — y K cos (n2); 


in the x direction it is — y K cos (nx). The surface integral in (2) 
accordingly becomes 


eer AMES oe, —xcos(nx)]dS. (58) 
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The integral here may be expressed as Z, + Z, — X; — Xz, where 
Z and X refer to the z and x terms in the integrand, and the subscripts 
1 and 2 refer to the two terms in K. Thus 


= SS = 2c0s(nz) dS, (59) 
1 
and similarly for the others. 
We find 
olds ra ™ cos? 6 sinddé 60) 
ORC hie i (a? — f? cos? 6) */? ’ ( 
which, with the substitution « = cos @ , becomes 
@203t0(2* oad 2 
Z, = 2a—— | ————_ 
1 fe if (Hay 97 (61) 
f? 


This may be evaluated from a standard table of integrals, e.g. Peirce, 
integral no. 157. 

Proceeding similarly, with the same substitution, Z,, X., and 
X, reduce to forms which can be evaluated from Peirce’s nos. 157, 
151, 139, and 127. The details are fairly lengthy and will be omitted 
here. The final result is 


Obs “ 
T=—azyb (62) 
at ai A 
— (3 —— 4—-) sin'-.- 
For a nearly spherical shape, ae which f/a is small, we find by 
expanding sin “ and b = a,| 1— — in powers of f/a, that the brack- 
2 
eted expression in (62) tends fnehe iB er . 
Thus j - 
5 gittong BRS Bix, ai 2% 63) 
6 Find Cid oid a8 esto bunt 
or, sinceea~ Dd, 
pe alappeiae =P (64) 
375 ab 


For very elongated shapes, with f ~ a, the bracket tends to a a/2 
and thus 
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ou oe tered nals (65) 
37 8 b 
Since’ a > > b this may be written as 
sh RATS yee) 66 
em a Oe Tat ee) 


If a = 2b the arc sine term in (62) drops out and the bracket is 
equal to a, which gives 


OPPS MA ta TG LO lard 3 2 (67) 


To consider whether é can change sign for intermediate eccen- 
tricities of the ellipsoid, we first note that the bracket is necessarily 
positive for a > 2b. It is zero for f = 0, and positive for f just 
greater than zero, as was seen in connection with (63). Hence (for 
simplicity taking a = 1) if it is to have a negative value at all it must 
vanish for some positive value of f less than \/3/4 (which corresponds 
toa < 2b). For this to occur there must be a root of the equation 


Ppt 68 
sun a5 V oS pam 2 ( ) 
between 0 and \/3/4. Now the derivatives of the expressions on the 
two sides are respectively 


1 
(left) VA ee (69) 
and 
(right) 9 — 24 f? + 32 ft — 16 f* 1 (70) 


9 — 24 f? 16 f4 V1—f? 

On the range 0 < f? < 3/4 expression (70) is positive and greater 
than (69). Hence the right side of (68) increases faster than the left 
and, since they are equal at f = 0, they cannot again become equal 
forife= iy) 3/4, 

Thus the bracket is always positive and the é contribution of the 
surface tension is always negative. This bears out the statement made 
above that a constant surface tension always opposes elongation. 


Rashevsky’s approximate expression for the surface tension ef- 
fect (1938, p. 311) is 


(71) 
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which with a: 6b = a: b becomes 


lda y §13a—bd 
fa ae <i) 
The above expressions may be summed up in the formula 
y a—b 
— ae a. eee (73) 
where the values of a are as follows: 
a: ellipsoid a: Rashevsky 
a=~b 1.6 1.15 
Cia 1.5 1.15 
eo 0 Lis iS 


The agreement is seen to be good, Rashevsky’s formula merely under- 
estimating the surface tension effect a little for nearly round shapes. 


Effect of variable surface tension: 


A surface tension which is not constant, but varies from point 
to point on the surface, is regarded by many physiologists (Wilson, 
192-197) as playing an important role in cell deformations. We shall 
here consider a simple case of this effect. 

For simplicity suppose that the variation is merely between poles 
and equator of the ellipsoid, and that it is given by 


» = constant + 7, cos’ 6. (74) 


Consider first the mechanical effect on the body due to the variation 
in perpendicular pressure on its surface. For this it is necessary only 
to amend the integral of (58) by inserting a cos? 6 factor into the in- 
tegrand. The procedure is then as before, except that a couple of new 
integral type forms are needed. These are given by Dwight (1934) 
as numbers 324.01 and 324.03. The final result for the contribution 
of this term to @ is 


OM LAI ot 
desmad ay.) Bins 2 
es, fe | 7 as f (75) 
Ls (eee 2 2 ry-1 i 
string ait Re ee a 
For a nearly spherical shape this becomes 
bt tenant Ge 8) 97 2. (76) 
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for a = about 24 0 it is 


gue le (77) 
n b 
while for a very elongated shape 
1 
5 oe a ee (78) 
yn 382 0 n b 
Thus in general 
——— ee 8 : ; (79) 
ie) 


where apparently £ increases from about + to 4 as the body elongates 
from spherical shape. A positive tendency to elongate requires that 
t, be negative; i.e. that y be greater at the equator than at the poles. 

However, the variation in normal pressure is not the only me- 
chanical feature involved. If the surface tension varies over the body 
there will also be tangential stresses present, since otherwise each 
surface element would be subjected to a net uncompensated pull. A 
simple illustration is a body covered by a sheet of rubber which is 
stretched non-uniformly. To maintain the non-uniform tension in the 
rubber it is necessary that friction with the body surface provide 
appropriate tangential stresses. But then the body is also subject to 
tangential surface stresses; i.e. the non-uniform surface tension gives 
rise to both normal and tangential stresses on the body. 

It might be possible, on the other hand, to have the requisite 
tangential stresses on the rubber supplied from the outside; in which 
case they would not be transmitted to the inner body. In general they 
may be partially sustained by the inside body and partly by some out- 
side agency. This latter would appear to be the case in actual cells, 
as illustrated by the experiments of Quincke (Wilson, p. 194) in 
which attendant streamings are observed both outside and inside the 
drop. 

For purposes of computation let us consider the first case, in 
which the entire unbalanced tangential surface stress is transmitted 
to the body inside. It is seen at once that the effect will oppose that of 
the varying perpendicular pressure, in that a y higher at the equator 
gives rise to surface pulls directed toward the equator and hence tend- 
ing to shorten the body. ; 

In the case of a system of revolution this tangential stress at each 
point of the body surface is directed toward the equator along the 
meridian curve and is equal in amount to d y/dl, where dl is the are 
length along the curve measured toward the equator. For an ellip- 
soid, and with y given by (74), this turns out to be 
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dy 2 t, cos 6 sin 6 
ee ee ee (80) 
dl VV a? — f? cos? 6 
Taking into account the angles between the surface tangent and the 
z and x axes we find the component surface stresses to be 


_ 2at, cos 6 sin? 6 


re a? — f? cos? 6 
(81) 
; 2 b t, cos? 6 sin 6 cos 
6 ea a IEEE LES? 
a? — f? cos? 6 
The surface element is 
dS = bsin6 \/a — f? cos? od6dd¢. (82) 


Upon writing out the expressions for ff z2Z'ds and ff « X' ds 
it is found that they reduce to one of the integrals encountered above, 
and that the final result is 


b 3 a? b 
ere yt ot a 
g= aot Aa (83) 
2 7) a b f a2 83 at J i 
teem ge ee sai: 
For a nearly spherical shape this reduces to 
Cig Wiha) 
ey 84 
ee (84) 
at G.— 2.b-itus 
Seeds One 
Se ect ak fom 85 
cee a he (85) 
and for long shapes 
t& 21 
6 ee ees : 86 
é soe (86) 
Thus the contribution of the tangential stresses is 
é= ids Oy 2 ? (87) 
ye oe 


where apparently a, is nearly equal to 0.4. 

Comparing with (79) shows the opposition of the two types of 
stresses as mentioned above. It should be noted that one or the other 
may predominate depending upon the shape of the body. For a nearly 
spherical shape the tangential stresses predominate so that a net aver- 
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age elongation requires that y be larger at the poles, while for longer 
shapes y must be larger at the equator. The quantities a and 6 have 
been evaluated above only at three shapes apiece, but from these it 
appears that the change in term predominance takes place at about 
a = 2b. This may shed some light on the controversy as to whether 
y should be larger at poles or at equator in order to divide a cell (Wil- 
son, p. 197). The presence of the constant term in (74) always acts 
to reduce é , as seen earlier, and so to oppose elongation. There is thus 
a variety of possibilities of elongated shapes having average length- 
wise equilibrium, in the sense that é is zero, under the action of a 
variable surface tension. 

If part of the tangential stress is sustained by the cell exterior, 
then (87) would be reduced and the normal pressure term (79) would 
become more dominant. 

The author is indebted to Mr. H. D. Landahl for checking and 
correcting portions of the calculations. 
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Using a recently developed approximation method, the problem 
of the deformation of a cell under the action of diffusion forces is 
studied for the general case when the rate of deformation is so rapid 
that the diffusion state is not quasistationary. The possibilities of 
configurations with several equilibria, as well as of periodical con- 
tractions and expansions of the cell around one of such configura- 
tions are shown. The bearing of these results on the theory of tissue 
contractility and amoeboid movements is discussed. 

In our book, “Mathematical Biophysics” (Rashevsky, 1938), here- 
inafter referred to as MB, we have given and discussed the equations 
of elongation of a dividing cell, under the influence of “metabolic dif- 
fusion forces”. We have considered in more detail the case when the 
elongation is so slow that at each moment the diffusion state may be 
considered as quasistationary. We have however given also, without 
discussion, equations for the general case, when the speed of deforma- 
tion may be arbitrary (MB, p. 315). The general equations have been 
corrected and amplified in a subsequent paper (Rashevsky, 1939). In 
the present paper we shall discuss somewhat more at length the case 
of deformations, which are sufficiently rapid so that the conditions of 
quasistationarity are not satisfied. In conclusion we shall point out 
the possible bearing of these cases on several important biological 
problems. 


rE 


Instead of equ. (45) p. 315 of MB, we consider the more general 
equation: 


re Cem aes. F 
Dera ee 
2(D. . 2a(aD. + 2bD;) -— 
—AnD;D- b?(De + 2Di) + 2a (aD. + 20Di) re = Co) A (1) 


(aD. 26D;) (De + 2Di) 


which is easily obtained by the same argument as equations (5) to 
(9) of a previous paper (Rashevsky, 1939), by putting h=od=b 
the notations here being the same as there. Putting 


aN ty ens 
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c a Co SG ’ (2) 
equ. (1) may be written: 
de 3D;D. b?(De + 2D;) + 2(aD. + 2bD;) a 
Gd ovy mabey. GD A=2bD) GaeseDs 


Instead of equation (47), p. 315 of MB, we have now the improved 
equation: 


(3) 


lda_ 2kTu 6D,D.(a — b) i 

adt  3My (26D;-+ aD.) (26D; -+ bD-) 
in which we have to put 6 = b (Rashevsky, 1989, equ. 23). Adding 
also to the right hand side of (4) the term —y(a — b)/3nab which 
expresses the effect of the surface tension (MB, p. 311, equ. 36), we 
obtain: 


(4) 


1da _ 2kT u D;,D.(a— b) ‘ y a—b (5) 
adt.-3My. (20D;4- aD.) (2Di4.De) 3n abs 


We shall now discuss the simultaneous differential equations (3) and 
(5). 


The requirement 
de 
abd eh | We 
7 =0 (6) 


introduced into (3), gives: 


_ ab’q = (aD. + 2bD;) (De + 2D) 
< 3D,D. 2a(aD.-- 2bDi) + O(Dsq 2Di) (7) 


The equality sign in (7) gives an equation of a line, for which dc/dt 
= 0. Since (MB, Appendix): 


3V 4 3V 
Chae Ta ; a te A = ie ; (8) 


the equation (7) may be written: 
7» -V4(De + 2Ds) Da + 2Di;V3V/4na 
4nDiD. 2D.a? +. 4D;\/3Va/4n-+ 3V (De + 2D) /Ana 


For very large values of a, we have: 


_ Vq(De+2Di) 1 
2a Sep (10) 


(9) 
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For very small values of a, we have: 


a 2q 3Va 
3D.N 42 


From (10) and (11) we see that c = 0 for a = 0 and for a = o. The 
value of ¢ is everywhere positive. From physical considerations it is 
clear that c should have a maximum for a = a = ¥/3V/4a, that is, 
when the cell is spherical (Fig. 1, full line). For D. = o this is ac- 
tually the case, as can be easily seen. For finite D. (9) gives a maxi- 
mum at a point slightly different from a,. This is due to the approxi- 
mations introduced by making 6 = b (Rashevsky, 1939) which is not 
quite exact. 
Because of equation (5), the requirement 


da 


(ir) 


ee 
Ti 0 (12) 
leads to 
_ My(2D,+ D.) 2D; | De 
Se ORT Dg nad sbi 3) 
for 
i> bi 
and to 
ae My (2D; — D.) 2D; D. 
cae TaD Deis eal LY >? uy 
for 


ab 
Because of (8), equations (13) and (14) may be written thus: 


My (2D; + ,) 2D; D; ao f b 15 

SRTDD. at a7 A) sit On pe 8 ieachy 1S} 
and 

_ My (2D: +D.) 2D; Daa : % RY: 

any sr mee a ar POLE a ae 

For a = b equ. (4) gives always da/dt = 0. The equality sign in 

(15) and (16) gives us an equation of a curve which is infinite for 

a = 0, varying for small values of a like 1/a ; while for larger values 


of a it varies as \/a (Fig. 1, broken line). It has a minimum for 
AD; 
Uy ere 2/8 q, . | (17) 


The two lines (9) and (15) may either not intersect at all (Fig. 
1) or they may intersect as in Fig. 2 or in Fig. 3. 
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FIGURE 2 
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The case of Fig. 1 will happen when for a given a, and other 
constants, the value of q is sufficiently small, resulting in a too low 


concentration ¢ = 


¢ — ¢ at the maximum, which corresponds to a 


sphere. In this case the two lines (9) and (15), together with the line 
AB which corresponds to a = a, divide the quadrant of the a,c — 
plane into six regions. In region I dc/dt < 0 because of (6) and (9), 
while da/dt > 0, because of (12) and (16). Any configuration point 
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a— 


FIGURE 3 


in the a,c — plane moves in the direction indicated by the arrow. In 
region IJ we have dc/dt < 0, da/dt < 0 ; the direction of the move- 
ment of a configuration point is also indicated by the arrow. In region 
III da/dt >0 ; de/dt < 0, while in region IV da/dt < 0 de/dt < 0. 
A similar survey of the other regions leads us to a distribution of 
pathways along which the configuration point moves according to 
equations (3) and (5) as shown on Fig. 1. For any point on AB, 
da/dt = 0. Therefore above the heavy line any such point will move 
downwards, below it — upwards. All paths converge to the point O, 
which corresponds to a = b = @y, that is, to a sphere. This is quite 
clear physically, since for very small values of q the sphere is a stable 
shape of equilibrium. 

A similar analysis of the case represented on Fig. 2 leads to a 
different distribution of the paths. The point O is now a point of un- 
stable equilibrium. There are now two points O, and O, of stable 
equilibrium, towards which the paths converge. The point O, corre- 
sponds to a < a, ora < 6 ; in other words to a flattened shape. The 
point O. corresponds to a > % or a@ > b, that is, to an elongated 
shape. This happens when, due to a sufficiently large q, the sphere 
becomes unstable for a given a), that is, when a becomes greater than 
the critical radius a* given by equ. (28) of a previous paper (Rash- 
evsky, 1939). In this case the metabolic forces begin to elongate the 
cell. This elongation may proceed either in one direction only, lead- 
ing to oblong shapes with an eventual division in two; or it may pro- 
ceed in two directions, leading to flattened shapes, like for instance 


red blood corpuscles. 
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From Fig. 2 we see, that while the paths of the configuration 
points converge towards O, and O;, they do not necessarily approach 
those points directly. These paths are in general convergent spirals, 
and represent therefore damped oscillations of ¢ and a around the 
points O, and O,. 

For the situation represented on Fig. 3 a possibility of undamped 
periodic oscillations exists. In this case both O, and O, correspond to 
a > a. But O, is unstable, while O and O, are stable. The paths 
around O, may be closed lines, representing undamped oscillations, 
which of course are not sinusoidal in general. If the case represented 
on Fig. (3) is actually realised, then a cell would be in a relatively 
stable equilibrium when spherical. It will require a finite external 
disturbance to elongate it to the length a, , corresponding to O,. But 
if a finite external disturbance brings the cell into a configuration for 
which a > a,, then the cell will stretch spontaneously further and 
will oscillate indefinitely, contracting and expanding rhythmically. 

The exact analytical determination of the paths of the configura- 
tion points presents a rather difficult problem, due to the nonlinearity 
of equations (3) and (4)> We may however study the behaviour of 
the cell in the immediate neighborhood of the points O, O, and O, , by 
expanding the functions v(a,c) = dc/dt and u(a,c)= da/dt in the 
vicinity of the corresponding points. But even then this requires an 
analytical expression for the coordinates a,c, and @2,¢. of O, and O., 
which it is difficult to obtain in an exact closed form. For a general 
orientation we may use however some rather crude approximations. 
As an illustration we shall use it for the study of the behaviour of the 
system in the vicinity of O, in the case represented on Fig. (3). 

That case will occur when the minimum of (15), given by 
(17), is shifted very much to the right. This happens when 
D; >.> D.. We shall therefore consider such a case, in which D, 
may be neglected as compared with D; , when the two occur as a sum 
or difference. Moreover we shall consider the case that a, is so far to 
the right, that instead of (9) we may use in this region the approxi- 
mate expression (10). The latter now becomes: 

fared 
~ AnD.’ feb, 


which is the approximate equation of the heavy line of Fig. (3) for 
large values of a. Subsequently we shall investigate the degree of 
approximation given by (18). ° 

Equ. (15) now becomes: 


My 2D; 


pe Preaea . 
Cie AD acta enaee ane ap aed, 
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The abscissa a, of the point O, is obtained by equating the right hand 
sides of (18) and (19), which gives: 

ee 1 My 2D; D. 
or rearranging: 


Vq oe MyD. (@ 
“6x RT us 2RT wu. de 


In MB, Appendix and a previous paper (Rashevsky, 1939) we have 
defined the critical radius of a cell as the radius above which the cell 
elongates spontaneously, for any deformation, no matter how small. 
We shall call the corresponding volume the critical volume V,. In the 
case represented on Fig. (3), V < V-, since between O and O, the 
spherical shape is the stable one, and the cell does not elongate spon- 
taneously at O. 


( =)22, (20) 


pey = a" (21) 


where a* is given by equ. (28) of a previous paper (Rashevsky, 1939), 
namely: 


9M (2D;+ D.)y 
SS sens ere i 
Combining (21) with (22) and neglecting D., we have 
12aM Diy 
= 2 
Vesa Vc RT aq (23) 
Putting 
ele! (24) 
Ve 
(20) and (23) give: 
(Gyan = art, 
or 
D; 
ya (ay , (25) 


While a may be negative, a. is always EEPORG because of the 2/3 


power. 
Introducing now the value (25) of a. into (18), we find the ordi- 
nate c, of O.. After some rearrangement this gives: 


Vaq 


= ———_—_. (26) 
Ana,%/D DD? 0? 
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Equation (5) reduces, fora > > 0, to 


da 2kT ue D,Dw eet i> (27) 
dt 3Mn (20D; + aD.) (2D; DA) we sab 


For D; > > D, and for very large values of a, this becomes 


da KTuD-. An 
pene alba oi —— g/2. (28) 
dt 6MnD; Vat 37 VN 3V 
Denote 
da dc 
=—__: =. 2. 
Lara ik = ae (29) 
Since, at O2 , u = O, we have in the vicinity of O2 
da ou Ou 
ee a Cr 0, (@ — 2) + ales (c— ¢), (30) 


where the subscript indicates that the values of the derivatives are to 
be taken at O., that is, fora = @&,c =. 


From (28) we have, using (8): 


Ou. 5RT uD. & yAe'/? 
25) GS 2 ga ee ee 
re % 12M7nD,a3" ses, 2703/2 (31) 
Because of (18), in which we must put now a = a@., c= ¢,, (31) may 
be written: 


Ou. =. pT AGV. ey ete. 
Gala ri (FaMD, ae dis! ; (32) 
Putting: 
5 V 
3V, 5 A 62 (33) 


(82) may be written, pecans of (23); 
OU. Gr 
(ae By 


oa ~~ Dna? (34) 


From (28) and (29) we also have: 
ES _ RTuD. 4n at”, 
oc’ 6MyD; BV" 
which, because of (8) and (25), may be written: 


OU, RT yd |a| ae ; 
oc 02 6My Gy —~) ee (35) 
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The absolute value |a| of a is taken in order to indicate that (35) 
is always positive because a, > 0. 
Introducing now (8) into (38), the latter may be written: 


dc _ 
apy q 
(36) 
4nD;,D-. 3V (D. + 2D;) /4na + 2D.a? + 4D;\V/3Va/4a 
oa (Qo 
CD, a 2D;) Da + 2D; V3V/4na 
For large values of a and for D; > > D, this becomes because of (21) 
De 
ove ac. (37) 
Since at O., v = 0, therefore in the neighborhood of O2, we have 
de ov ov 
a ae ae (a — a2) + (0s (¢ — C2). (38) 
Using (8) and (25) we obtain from (37): 
ov Ga; A ov 3~0/ DD? 0? 
eee ee (ee en (Ran we eee ee 39 
(Gor Qo Cp ? Co) 0s a, ( ) 
Introducing now the notations 
e=A—G; Y=C—G; 
ByQ,'/? = +94) ( Ki ) 23 ; 
ero ae @y) | aD; : 
(40) 
F 3, 2.72 
ee sTe Ee Pye nie 
6M D. a" 
we find from equations (30), (88), (29), (34), (35), and (39): 
dx 
a AS B, 
Ti a+ By 
(41) 
dy 
—— = Asr--B; 
Ti 2t + Boy 


The solution of this system of equations is well known (Rashevsky, 
1938a). It possesses undamped periodical solutions, if 


A, + B,=0 (42) 
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and 


A,B, —B,A,>0. (43) 


Requirement (42) together with (40) gives, after simple rearrange- 
ments, using (25): 


6D22D in es B 


: (44) 


YU Vy |a| 


From (24) and (23) it follows that in general B © a~©1. Tak- 
ing as plausible values D. ~ 10° cm® sec* ; Di ~ 107% cm® sec” ; 
a) © 10-? em., (44) requires that n/y ~ 10-° cm™ sec., which is in 
agreement with a result found elsewhere from different considerations 
(Rashevsky, 1938c). 

Requirement (43) together with (40) gives after rearrange- 
ments, using again (25): 


96MD; 
a, ae (45) 
RT uq 
which on account of (22) may be written (remembering D; > > D-) 
ao? > Bats. (46) 
Inequality (46) requires: 
a 3 
> B. (47) 
But, because of (33): 
Ao° Vee 2 
ge ay ae ED: (48) 


Hence (46) is equivalent to: 


2(6+1) >28, 


or 
24 
| B< 3° (49) 
Inequality (48) gives now: 
iY $792 
V. =< 3? (50) 


and this, because of (24), makes 


a<0Q. (51) 
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Periodical contractions and expansions may thus exist under 
rather general conditions. Taking V/V, = 1/3, which makes a = 
—1, and which fixes qg because of (23), since V is fixed by the choice 
of a , we find that a. = 465 a. 

Putting 

gan; D, = nD.» (52) 


Qo 
the second fraction in the right-hand side of (9) may be written 
1 22+ Qnv/z 
Qy 222 + 4nzV/z+ 2n-+1 

In our case we have n = 10 ;z = 4-65. The first term in the de- 
nominator of (53) is equal to 200, the second to 400 and the last two 
to 21. The omission of the last two terms introduces an error of about 
A%. Omitting (2n + 1), (53) may be shortened by z? -+ 2n\/z and 
becomes equal to 1/2a,.z = 1/2a. 

This leads to equ. (10) as an approximation to (9), and we see 
that, with the values of the constants used here, the approximation 
is good. 

With the values of the constants used above we find for the fre- 
quency 


(53) 


w= VA,B,— B,A2 


an expression inversely proportional to the square root of 7. With 
n . © 10° gmcm™ sec” . which is a very high value, corresponding 
to a solid jelly, we have w © 10° sec”, giving a period of the order of 
hours. With smaller values of 7 the period can be made much smaller. 
With slightly different values of the constants we may have: 


A,t+tB,>0; (A, + Bi)? — (A,B, — B,A2) < 0 (54) 


=—_———_. 
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FIGURE 4 
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which will give oscillations with an increasing amplitude. The ampli- 
tude of course cannot increase indefinitely. What happens 1s shown on 
Fig. (4). In the neighborhood of O, the amplitude increases (broken 
line). Further away it decreases (broken line). Somewhere in be- 
tween lies the stable closed path (full line), corresponding to un- 
damped non-sinusoidal oscillations. An exact analytical treatment 
of the case presents a problem of great interest but of much greater 
difficulty than the present one. 

It is of great interest to study possible cases in which the length 
of the cell oscillates periodically around an equilibrium value, so that 
the cell periodically contracts and expands. Such studies may even- 
tually form a basis for the theory of amoeboid movements, or periodi- 
cal contractions of tissues (heart). 


II. 


In Chapt. VI of MB we have shown on a limiting case that, in 
case of two interacting substances in a cell, periodical variations of 
their concentrations are possible. A. Weinberg (1938) has given a 
treatment of the problem for a more general case, using the exact 
solution of the diffusion equations for a sphere. Later on (Weinberg, 
1939) he obtained the same results by using our new method, appli- 
cable to non-spherical shapes, developed in MB, Appendix. 

If the concentrations and their gradients vary periodically, so 
will the diffusion forces vary also, and it is natural to expect that un- 
der the action of such periodic forces the deformation of the cell also 
varies periodically. A method somewhat similar to the one used above 
may be employed for treating that problem. We shall now have two 
substances, the average concentrations of which we shall denote by 
c’ and c”. The external concentration will be ¢,’ and ¢,”. and all 
other constants referring to the two substances shall be marked by 
primes and double primes. We do not use subscripts , and » like in 
Chapt. VI of MB, because those subscripts have a different meaning 
in the derivation of equation (1) as given in a previous paper (Rash- 
evsky, 1939). Using a method similar to Weinberg’s (1939), we can 
now determine, for a prescribed a and b, the values c”® and c”, 
around which c’ and c” will oscillate. These values, plotted against a, 
will give us curves that are zero for a = 0 and fora = ow. Due to the 
additivity of the diffusion forces, an equation, corresponding to equ. 
(5) will give us da/dt as a linear combination of ¢’ and c” 


da ; 
qe Ae + Be" +C, (55) 


where A, B and C are functions of the reaction rates, diffusion co- 
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efficients, surface tension, etc. The requirement da/dt = 0 gives us 
then Ac’ + Bc"as a function of a: Ac’ + Bc” = F(a). Since we also 
have ec’ and c”® as functions of a: c’° = F(a); c"° = F,(a) ; from the 
Fequirement dc’'/dt = 0 and de” /dt = 0, we can plot Ac® + Bc”? = 
AF, (a) + BF,(a@) on the same diagram. The root of F(a) = AF,(a) 
-+ BF.,.(a) gives then the values a, of a, for which de'/dt = de” /dt = 
da/dt = 0. The corresponding values F,(a,) and F.(a.) give us the 
values of c’, and c”,, corresponding to that configuration. We may 
then investigate the behaviour of the system in the neighborhood of 
that configuration. 


Lif. 


The assumptions made in MB, Chapt. VI concerning the inter- 
action of the two substances, which lead to equ. (20) and (21) are 
somewhat artificial. From the point of view of chemical kinetics it 
seems to be more plausible to assume, that the rate of production of 
the first substance is given by 


Oy (A, = ic”) C ’ 
indicating that the second substance inhibits the formation of the first. 
Similarly we may take for the rate of reaction of the second substance 
— a, (Bs — y2e')e", 
a,, 61, y1 etc. are constants, which are not to be confused with o., B 


and y of section I. 
Instead of equ. (3) we now have: 


ie = a,(6:— pie”) 

(56) 
3D,'D. b?(D- + 2D;') +2 (aD, + 2bD;') (c ra 
ab? (aD, + 2bD,') (D. + 2D’) 

ae = —o2 (B2 — y2€') c" 

(57) 

3D;"D." b?(De-" + 2D;”) + 2(aD.” + 2bD;") ee 


man gb? (aD DOD 2) (Dre 2Dp yy 


- Consider the special case, when ¢,’ = ¢)” = 0 and when the cell is im- 
permeable to the second substance, so that the latter does not flow 
either outward or inward. For instance D,” = 0 ; then 
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de’ 
Gt = [oi6, 
(58) 
SD/Dé_b4(Dé + 2Ds') +2(0Dé + 20D) aa, 
——— ee TD Tv Ne ee ? 
ab? (aD, + 2bD;’) (De + 2D;') 
ca = — a2 ($2, — yo’) ¢ ; (59) 
dt 
which for a given constant a are of the form 
ue = 
a = ¢' (A, — pe”) 5W 
(60) 
dc” So. =) 
‘7 ies Oh Ase 


Vito Volterra (1931) has shown that equ. (60) have nondamped peri- 
odical nonsinusoidal solutions under rather general conditions. It is 
therefore interesting to investigate whether the more general system 
(56) and (57) has similar solutions and then pass to the still more 
general case, by adding to the system (56) and (57) an equation 
corresponding to equ. (5). 

The great biological interest of this type of equation lies in the 
fact that, whereas the linear cases have periodical undamped solu- 
tions only for rather restricted values of the constants, equation (60) 
have them always for positive values of 24,, 42, uw, and ws. Moreover 
the solutions, while periodic, are not simple sinusoids and represent 
therefore a much more general type of oscillations. Volterra (1931, 
p. 56) has studied still more general systems, which have undamped 
oscillating, but not periodical solutions. 


Vic 


The considerations of this paper bear on two different but re- 
lated biological problems. We have here a somewhat more general 
theory of nonspherical cells, than the one given previously (Rash- 
evsky, 1938b). In the latter we required the vanishing of all princi- 
pal stresses, a requirement which is not necessary for stable equilib- 
rium. See also the correction to that paper, mentioned later (Rash- 
evsky, 1939). Going back to section I, we remark that the equilibrium 
configuration O, is a function of qg, D, T and y. In the general case 
it will be also a function of h. A change of any of those parameters 
results in a change of a, , hence in either a contraction or elongation 
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of the cell. Such a contraction or elongation may therefore ‘result 
either through a change of metabolic activities of the cell (change of 
q), or through a change of permeability (D and h), or through both. 
The same will hold for the more general cases, discussed in sections 
II and III. Moreover, things will be complicated by the possibility of 
periodic contractions and relaxations, which under circumstances may 
go on indefinitely. The importance of the above for a possible physico- 
mathematical theory of contractile tissues is evident. We do not wish 
to imply that the simple model discussed here is even remotely 
adequate to represent a real muscle-cell. But, in all our previous 
studies, we must begin with rather abstract oversimplified cases, 
which may not apply immediately to actual cases (MB, introduction). 
From this point of view the above considerations may be considered 
as laying down the foundations for a biophysical theory of contrac- 
tility, which incidentally has the advantage of following naturally 
from our general mathematico-biophysical theory of the cell, in the 
course of a systematic development of the latter. 

The actual force developed during the contraction of a cell studied 
here can be estimated by evaluating for instance o,. From equations 
(14) and (17) of MB, Appendix, we have, for D. = w: 


_2RTyq ab? 
- 3 MD 2a?-+ 0 
with the values of the constants used in MB, and witha ~ b ~ 10° 
cm, we have: 


(61) 


oC 


o, © 10? dyn-cm”. 


A threefold increase of a, and.a tenfold decrease of D , which leaves 
the latter well within plausible limits (MB, p. 53). would give: 


a, © 10* dyn-cm? = 10 gr-cm” . 


A finite D. does complicate the expressions, but does not change the 
order of magnitude. 

Another biological phenomenon on which the foregoing considera- 
tions may have a bearing and the possible connection of which with 
contractility of tissues has been already pointed out (O. Furth, 1922) 
is that of amoeboidal movements. Our present discussion gives us a 
picture of a cell, which expands and contracts in one direction, not 
necessarily periodically. If noticeable inhomogeneities are present in 
the cell, these will produce strong asymmetries in the distribution of 
the forces, and may result in elongation in several directions and for- 
mation of protrusions like pseudopodia. The study of the above men- 
tioned types of equations, which give oscillating but not periodical 
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solutions (Volterra, p. 56) should be of particular interest in this con- 
nection. It is also to be remarked, that any nonuniformities of the cell, 
due to coarse inclusions, will be subject to statistical fluctuations, 
which will impart a degree of irregularity to the movement of the cell, 
If we have N coarse particles per element of volume, the average 


fluctuation from the mean is equal to 1/\/N (Gyemant, 1925). With 
1000 coarse particles distributed at random in a cell, the average fluc- 
tuation in a volume equal to 1/10 of the cell volume, and containing 
therefore on the average 100 particles, will be 10%. This will affect 
the distribution of the gradients and forces by about as much and will 
result in rather irregular movements. This leads to an interesting 
problem of the combined effects of regular periodicities and of sta- 
tistical fluctuations. 

The author is indebted to Dr. A. S. Householder for checking the 
calculations. 
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CONCERNING RASHEVSKY’S THEORY OF THE “GESTALT” 
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In experiments of Kliiver on monkeys and of Hertz on bees, 
observed reactions indicate the presence of a neural mechanism which 
serves to order geometric figures in a linear sequence, so that any 
given figure occupies a definite place in a one-dimensional array. In 
a theoretical mechanism described by Rashevsky, presentation of a 
visual Stimulus-figure results in an excitation at a given center whose 
intensity is a function of the contour of the figure. This functional 
dependence of excitation upon the curve is examined in some special 
cases, and the mode of variation is further described qualitatively. 


I. The Problem 


In the course of an extensive series of experiments performed by 
Kliiver (1), a monkey is presented with two boxes, identical except 
that the nearest face of one has on it a large rectangle, while that of 
the other has a (nearly) similar rectangle of half the area. By put- 
ting food in the first box and none in the second during a period of 
training, the monkey is trained to react “positively” to the larger rec- 
tangle. After a period of such training the pair of rectangles on the 
boxes is replaced by another pair having the same relative but dif- 
ferent absolute sizes. It is found that, unless the absolute sizes are 
made very small or very large, the monkey continues to react posi- 
tively to the larger rectangle, even though this may be identical with 
the smaller of the pair used in training. Moreover, if the pair of rec- 
tangles is replaced by a pair of figures having any one of a variety of 

different shapes, similar to one another and with areas in the ratio 
approximately 2 to 1, then the monkey still chooses the larger figure. 

This type of experiment is varied in a number of different ways 
and the results do not seem intuitively surprising. However, a some- 
what more mystifying outcome is obtained when pairs of different 
shapes of approximately the same area are utilized, e.g., a circle and 
a square. After the training with this pair of stimuli, the circle and 
the square are replaced by some other pair of distinct shapes. In gen- 
eral one of these shapes was chosen to give a certain impression of 
roundedness, the other of angularity. Then through quite a wide 
range of such stimulus-pairs the number of positive reactions to one 
was clearly in excess of the number of positive reactions to the other. 
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On the other hand, with certain other pairs the reactions were equally 
divided. b 

In the first experiment mentioned there is present with each 
stimulus-pair an obvious ordering relation, the same for all, with re- 
spect to which the two stimuli in each pair can be ordered. This rela- 
tion is evidently that of being-greater-than, and a priori it might 
equaily well refer to areas or to perimeters. We may therefore think 
of the monkey’s training as conditioning him to react to the first mem- 
ber of the ordered pair. A theoretical mechanism capable of this type 
of conditioning has already been described (3, ch. XXVI). 

The results of the other experiment can be made intelligible, when 
viewed in the light of the first, provided we can exhibit a common 
ordering relation of such a nature that the order is clear in the pairs 
between which the monkey’s choice was unequivocal, but obscure 
where the choice was random. Concretely, we require an arrange- 
ment of neurones, terminating in a single neurone or set of neurones, 
and transmitting to each terminal neurone an excitation whose inten- 
sity is a function only of the shape of the visual image. Then if the 
initial training conditions the animal to react positively to the shape 
yielding the excitation, say, of greatest intensity, we should expect for 
any other pair of shapes a clear preference again for the shape of 
greatest intensity, if these intensities are significantly different, but 
a random choice when they are (nearly) the same. 

Experiments performed by Hertz (4) on bees exhibit a fairly 
well-defined order of natural preference among the numerous shapes 
presented, quite apart from experimental conditioning. In general the 
bees prefer an abundance of detail within the total area covered, and 
they prefer small figures to large. However, Hertz distinguishes be- 
tween “figural intensity,” dependent upon the total perimeter, and 
“figural quality,” and finds a preference often clearly exhibited for a 
shape of one figural quality over another of different quality but equal 
intensity. Thus the intensity of excitation, on our theory, is a func- 
tion both of Hertz’s figural intensity and of figural quality. 

Rashevsky has demonstrated the possibility of such an associa- 
tion between shapes and intensity of excitation by his fruitful picture 
of a succession of intra-inhibiting synaptic levels (3, esp. ch. XXVII, 
XXVIII). The excitation is relayed from the initial sensory to the 
ultimate motor level, the transmission from any level to the next be- 
ing effected by a group of neurones whose arrangement determines 
the particular features of the distribution that are to be kept invari- 
ant. 

Now consider any one of the stimulus shapes S. Suppose the 
illumination, color, etc., to be uniform and the same for all, and as- 
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sume, in fact, that only the contour C, is illuminated. With a suitable 
arrangement of preceding neurones, the excited synapses at some level 
will be arranged in a curve which is geometrically similar to the pre- 
sented contour C,. This is, therefore, a closed curve C. Suppose the 
synapses comprising this level are uniformly distributed and numer- 
ous enough to be regarded, with sufficient approximation, as mathe- 
matically dense. 

Now suppose that each excited synapse inhibits each other syn- 
apse proportionally to its own excitation and to a function of its dis- 
tance from this other synapse; then if s and ¢ are the arc lengths 
measured along C from some fixed point of C to two given points, if 
r(s,t) = r(t, 8) is the ordinary distance between these points, and 
if 6(t) is the resultant excitation per unit arc length at ¢, then the 
inhibiting effect of t upon s can be written 2 K(7)¢(t), where J is 
some constant and K some positive function of 7. Then the total in- 
hibition acting on s is 


ASK (r)¢(t) dt, 


where the integration is made around the entire curve C. Thus the 
resultant steady state excitation (per unit length) per unit of im- 
pressed excitation is given by the solution ¢(s) of the integral equa- 
tion (3, ch. XXVII) 


é(s) =1—AsK(r)¢(t)dt. (1) 


If all synapses at this level are connected to a single synapse at 
the next level, then the total excitation at this synapse will be meas- 
ured by 


=> f4(s)ds. (2) 
For a fixed inhibitory function K , @ is a functional, dependent upon 
the particular contour C about which excitation takes place. The re- 


lation between the values of ® for contours of different shapes fur- 
nishes the ordering relation sought. 


II. Some Numerical Approximations 


We wish to investigate the dependence of the total excitation ® 
upon the shape of the curve C,, when the inhibition function K is a 
given function of r. Evidently the curves C should be taken with con- 
stant length. Actually, in any experimental check, the area, too, 
should be kept fixed. In the present preliminary investigation, how- 
~ ever, we keep only the perimeter constant. It is no restriction to take 
this perimeter as the unit of length, so that the limits 0 and 1 should 
be put on the integrals in (1) and (2). 
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We note that the simplest curve of all is the circle, for which 
the excitation is always constant whatever function K (7) is used. The 
solution of the integral equation can then always be obtained in closed 


form. Let 
ae [K(ryat. (3) 


Then k is independent of s for the circle, and the solution to (1) is 
¢=(1+41k)*. (4) 


Since the circumference is of unit length this is also the value of @. 
For a fixed &, ® vanishes asymptotically as 4 increases, and for all 
positive A, ® is less than unity. 

A priori it seems plausible that the inhibiting effect should be 
greatest between nearest centers, and that K , therefore, should be a 
decreasing function of r. A single increasing function was tried in 
two cases, viz., K = r?, for the circle, and for a double line-segment 
regarded as a rectangle of zero area. For these two cases, with 4 = 1, 
® turned out to be .9518 and .9603, respectively, and this small dif- 
ference for such different shapes was regarded as further evidence 
that an increasing function was uninteresting for this investigation. 

A natural choice for a positive decreasing function is r?. How- 
ever, checking with the circle shows that for a power of 7 whose ex- 
ponent is —1 or less, the solution does not exist. Evidently the in- 
hibiting effect could not become actually infinite. But assuming the 
effect to vary approximately as rr, though never exceeding some 
finite value, the function finally chosen was 


Ke (eee) oe 1 (5) 


The constant a was taken as of the order of magnitude of the distance 
between adjacent synapses. 

As a closed solution cannot be hoped for in general, a method of 
approximate solution was employed similar to a method described 
by Nystrom (2). The Neumann series solution of the integral equa- 
tion giving the solution as a power series in J is obtained by first 
replacing the ¢ in the integrand by ¢1 = 1 (the impressed excita- 
tion in our case, which is constant) to get ¢é2, then replacing ¢ by 
$2 to get ¢;, and so on. But these quadratures themselves become ex- 
cessively complicated if not impossible, and so are approximated by a 
suitable method of numerical quadrature. Thus if one tabulates the 
values of K for a sufficient number of (s , t)-values, and for each s 
calculates the quantities K(s , t)¢:(t), then the numerical quadrature 
gives the approximate value of $i.:(s). 
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If, for each s and7, K(s, t)¢(t) were continuous and had a con- 
tinuous derivative in t, then the most suitable method of numerical 
quadrature would be Gauss’s or Tchebyscheff’s. Actually, however, 
for 0 < s < 1 this function of f rises to a single maximum, on each 
side of which the curve is concave upward, and at which the deriva- 
tive is discontinuous. Therefore Simpson’s rule seemed to be the only 
feasible method. But again Simpson’s rule approximates to the curve 
in question by a set of parabolic ares, each are being determined by 
three consecutive points of subdivision. Evidently if the peak point 
were the middle point of a triple determining such an arc, this arc 
would give a very poor approximation. Hence, whenever the peak 
point entered as an even point of subdivision (counting the end-point 
as the first) the first and last arcs were approximated by line-seg- 
ments, and Simpson’s rule was employed on the intermediate portion. 

Further consideration made it possible to lighten the labor some- 
what. Evidently ¢: = 1 is not the only possible choice of function with 
which to start. In fact, suppose that $, and ¢, are two functions satis- 
fying 


11> >h>¢ 
for alls. Then, since K > 0, we have 
O<&ihe<e, 


where ¢. and $2 are obtained by substituting $, and ¢, respectively in 
place of ¢ in the right member of (1). Moreover, if 


1>$%>¢>H%>0, 
gir hrt, 
where ¢,, ¢, and ¢, are otherwise arbitrary, then 
oo Soe < $e» 


Now ¢: = 1 is everywhere too large, and the resulting ¢. therefore 
too small. But if any approximation were known in advance to be 
better than ¢. = 1, then the corresponding sequence of approxima- 
tions would stabilize more rapidly. In practice it appeared advan- 
tageous to start with a suitable ¢:, obtain ¢. by quadrature, and take 
for ¢3; some weighted mean. ' 

Consider now some actual results, starting with the circle for 
which the solution is exact. We find from (3) and (5) that 
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kes A LBD e 
so that (4) gives 
¢=6= (14 4.1835 1)>. (6) 
For a few values of 1 we have 
eae O:05% 0025 0.005 
(7) 
& = 0,193, 0.827 , 0.923 , 0.980 


Note that the successive steps of the numerical method outlined 
above, if applied here with ¢, = 1, correspond to the sums for increas- 
ing n of the first n terms of the expansion of (6) in powers of 4. How- 
ever, this does not converge except for 4 = .2 or less. The significance 
of this becomes clear when we consider for a moment the physical 
picture. When excitation first occurs it is momentarily uniform and 
free of mutual inhibition. As nearby centers have time to act, if the 
inhibiting effect is too great there will result a complete annihilation 
of excitation. But this would also destroy the inhibition and allow 
a re-establishment of excitation. Such an oscillation might go on in- 
definitely without a steady state being possible. 

To follow the dynamics in complete detail, would require some 
assumptions, or information, about the time-rate of inhibiting action. 
Hence we shall assumed to be small enough to permit the rapid estab- 
lishment of a steady state. In fact, we shall assume it small enough 
so that the initial inhibitory action nowhere completely inhibits. Hence 
every approximation ¢ > 0 everywhere on the curve. 

The next curve considered was the double line-segment, regarded 
as the limit of a rectangle of dimensions (144 — «) & ¢. From sym- 
metry, 


o(s) =9(1/2+s) =¢(1—s). (8) 


Hence we need consider only the interval 0 < s < 4, and the equa- 
tion can be transformed to 


#() =1—24 [ "(518474 (01 —|s—t)) Na(Hat. (9) 


The first quadrature required in the Neumann series can be carried 
out exactly to give the first approximation 


é:1(s) = 1— 22 log(1-+ 100s) (51 —100s) . (10) 
Ten subintervals were used. 
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For 4 = .005, we find 


t Py Py 
.000 .9607 .9614 
.025 9487 .9448 
.050 94388 .9445 
.075 9409 .9376 
.100 .9889 .9401 
125 .93875 .93847 
-150 .9365 9379 
175 .93857 .9830 
.200 .9852 .93869 
22) .9849 .9838 
250 .93848 .9365 
1 1 

{f gidt = .9399 , f godt = 9884. (11) 
0 0 


As we have seen, ¢, must be too small and ¢2, if computed ex- 
actly, would be too large. Hence our values for ¢2 cannot be accepted. 
The true value of & is somewhat in excess of .940, and better results 
could be obtained only by increasing the number of subdivisions. 

For 2 = .02, several successive substitutions and means were 
taken. The last two sets obtained are: 


t on Pnsa 
.000 -8675 86738 
025 8193 8177 
.050 .8106 .8105 
.075 .8018 .8009 
100 -7963 -7966 
125 -7922 7925 
.150 -7894 -7900 
175 -7872 7877 
.200 -7860 .7869 
225 7789 .7786 
250 -7850 .7861 


f “bn dt = .7978 , if "“buea dt = 7977. (12) 
0 0 
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For 14 = .05, we found 


t $5 ?, 
.000 Tol .7418 
.025 637 .6491 
.050 .622 .6386 
.075 .608 .6240 
.100 099 .6162 
125 .593 .6099 
.150 .588 .6062 
175 .585 .6030 
.200 .583 .6015 
229 581 .6002 
.250 581 .6001 
1 1 
f bs dt = 604, { br dt = 6204. (13) 
0 0 


Here ¢, was an arithmetic mean, and ¢, obtained from the quadrature. 

For 2 = .1, ¢, is negative except near the ends of the segment. 
The calculations were therefore not continued. 

As with the circle, @ decreases rapidly as 2 increases. For a fixed 
A, ¢ increases rapidly on approaching the ends of the segment. On 
approaching the midpoint of the segment ¢ decreases, but very slowly. 
This is pretty much what would be expected intuitively. 

The next curve to be considered was the Greek cross with limbs 
of negligible width. Since for this 


¢(s) =¢(1/4+s) =4(1/2+38) 
= ¢(38/4+ 8) = ¢(1— 8), (14) 
we consider only the interval 0 < s < 1/8, and reduce the integral 
equation to 


6(s) =1—24 f'"(Lo1+ |t—si}>4 [26—t—s} 


+ 2[.01+ ((4—s)?-+ (4—4)}']4} o(t)dt. (15) 


For 4 = .05 the first quadrature gives negative values of ¢2. We con- 
sider only 4 = .02. 

The method employed here was altered slightly by starting with 
only a single subdivision and then adding more as a fair degree of 
stability was obtained. The single subdivision gave ultimately: 
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t te e; 

0 .8308 8295 
1/16 -7015 .6983 
1/8 .6418 .6365 


Interpolation gave approximate values of ¢, for t = 1/32 and 
t rd 3/32 , and the process was carried on with four intervals of sub- 
division. Ultimately we found: 


t P10 Pu 
0 8465 8417 
1/32 -T785 -7765 
1/16 -7626 -7602 
3/32 -7252 £7222 
1/8 .6782 6742 


A final subdivision gave: 


t $3 Si4 

0 8457 .8465 

1/64 .8083 .8095 

1/32 -7894 .7953 

3/64 -7770 .78388 

1/16 7671 .7685 

5/64 -7518 O47 

8/82 -7345 -7389 

7/64 11138 rats | 

1/8 .6859 6895 
E 1 

f dus dt = .7628, i daadt-= 57664: (16) 

JO 0 


Noting that increasing the number of subdivisions increases the val- 
ues of ¢ at stabilization, we may regard ¢., as fairly accurate. 
These results, rough and incomplete as they are, furnish neverthe- 
less some indications as to the type of effects one may expect from 
such mutual inhibitions. Comparing the three contours, circle, rec- 
tangle of zero area, and Greek cross of zero area, for the case 4 = .02, 
we see that the constant excitation around the circle is greater than 
the maximum excitation for either of the other two. At the end of the 
rectangle and at the end of the arm of the cross, the excitations are 
not greatly different, whereas at the middle the excitations differ by 
about one-tenth. The total excitations for these two contours are about 
80 and .77 respectively, as compared with the value .92 for the circle. 
One must observe at this point that these total excitations are applied 
to neurones whose thresholds may be not greatly different. Hence ef- 
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fective differences between these excitations are not within a range of 
12 per cent or 15 per cent as they appear at first glance, but may well 
run up to several hundred per cent when, say, the neurone has a thres- 
hold of .75. : 
In general we may expect “corner” and “overlap effects of di- 
minished excitation near corners and points of high curvature, and 
also on arcs which lie close together. Hence decreasing the area, and 
so bringing the arcs closer together, should tend, in general, to cut 
down the total excitation. And again, cutting down the maximum dis- 
tance between pairs of points, even though the area remains the same, 
amounting to increasing the curvature and the number of corners, 
would diminish the excitation. Thus “jagged” curves, with numer- 
ous corners and crinkles, will yield low excitation, while rounded, full 
curves, ellipses with small eccentricities, will have a higher excitation. 


III. Some General Considerations 


An analysis in general along some such lines as the following 
seems indicated. Let K(x) be a properly monotonically decreasing 
function of r possessing suitable continuity properties. Consider the 
integral equation (1) associated with any closed curve C of unit 
length, which is continuous and whose derivatives have at most a finite 
number of points of discontinuity. If C does not intersect itself, then 
it bounds a region of the plane of a certain area. We may thus in- 
quire, 

(a) What curve C enclosing a given area gives to ® as defined 
by (2) the largest possible value? One would expect.to find a unique 
curve C,, bounding a region of given area, A, for which 6 = 4, isa 
maximum. If so, we may inquire 

(b) Is 4 an increasing function of A? It seems likely that C, 
would turn out to be everywhere convex, possibly an ellipse, and that 
6, is increasing. If so, the circle would furnish an absolute maximum 
for all A, since it bounds the greatest possible area. For a given 
curve C which encloses an area A and is different from C,, one may 
then hope to find all curves enclosing the same area and having the 
same total excitation 6 as C. 

Similarly one could ask for a curve with given area giving 6a 
minimum value. 


In making such an investigation it would be necessary to require 
that 


se cg 


in order to insure that a steady state can be reached for any curve 


A. S. HOUSEHOLDER 713 


whatsoever. With this condition the sum of the first two terms of the 
Neumann series will be positive at every point of any curve. 

It is hopeless at present to attempt any detailed check of this the- 
ory against experimental results. It has been the aim of this presen- 
tation only to show that some linear ordering is possible, and to indi- 
cate for a particular method of ordering [a particular K (17) ] how the 
figures would be ordered. However, we may note that in the experi- 
ments of both Kliiver and of Hertz, the contrasting features are 
smoothness as against angularity, simplicity as against detail, and 
this contrast precisely is obtained by the sharply decreasing functions 
K(r) here used as a measure of the inhibitions. Furthermore, it is 
just such a function which seems indicated from simple physiological 
considerations — nearby centers should affect one another more 
strongly than distant ones. Replacing the kernel (a -+ 7)-* by an- 
other function K(7) should alter only details, not the general char- 
acter of the ordering. Hertz’ results, then, seem to indicate a prefer- 
ence for figures giving lower excitations. 

Evidently there is involved not only the pure form, but also the 
size. The important element in the size may be the total perimeter ; 
it may be the area enclosed by the contour; it may be the area of the 
smallest convex region enclosing the contour; it may be a combination 
of all of these. Quite likely the ratio of these two areas, as a measure 
of the complexity, would be important. At any rate, any experimental 
checking of such a mechanism requires keeping constant as many of 
these elements as possible while varying the form. 
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